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Abstract 

The Reeb graph is a construction which originated in Morse theory to study a real 
valued function defined on a topological space. More recently, it has been used in various 
applications to study noisy data which creates a desire to define a measure of similarity 
between these structures. Here, we exploit the fact that the category of Reeb graphs is 
equivalent to the category of a particular class of cosheaf. Using this equivalency, we can 
define an ‘interleaving’ distance between Reeb graphs which is stable under the perturbation 
of a function. Along the way, we obtain a natural construction for smoothing a Reeb graph 
to reduce its topological complexity. The smoothed Reeb graph can be constructed in 
polynomial time. 
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1 Introduction 

1.1 Purpose 

The Reeb graph, originally defined in the context of Morse theory [35], can be used to study 
properties of a space through the lens of a real-valued function by providing a way to track 
and visualize the connected components of the space at levelsets of the function (Figure 1). 
When an algorithm for computation was given in [36] , the rediscovery of the Reeb graph by the 
computer graphics community immediately showed the Reeb graph to be an extremely useful 
tool in many applications. These include shape comparison [23,28], data skeletonization [12,25], 
surface denoising [47], as well as choosing generators for homology classes [18]; see [6] for a 
survey. Two main properties of this construction have made it extremely useful in the applied 
setting. First, its dependence on the chosen function and not just on the space itself means 
different functions can be used to highlight different properties of the underlying space. Second, 
it is rather quick to compute and thus can be used for very large data sets [20,27,32], 

Much of the literature is dedicated to studying Reeb graphs in the context of Morse functions 
where a great deal can be said about its properties [2,14,15]. In addition, several variations 
on the Reeb graph have also been proposed and proven quite useful. One such is Mapper [37], 
which applies the ideas of partial clustering to Reeb graphs in order to make the construction 
more robust to noise; this has found a great deal of success on big data sets [31,49]. A similar 
variation called the a-Reeb graph was used in [12] to study data sets with 1-dimensional structure. 
Another variation is the Extended Reeb graph [7, 23] , which generalizes the theory to non-Morse 
functions. Finally, a Reeb graph defined for a function on a simply-connected space cannot have 
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any loops, and is called a contour tree [10,33,41,45]. Any such contour tree can be replicated as 
the contour tree of a function on a 2-dimensional surface, thus allowing for informed exploration 
of otherwise hard-to-visualize high-dimensional data [26,46]. 

Many of these applications involve data, where one should operate with at least a modicum 
of statistical integrity. Therefore it is important to consider not just individual Reeb graphs, but 
the whole ‘space’ of Reeb graphs. In this paper we will: 

• define a distance function between pairs of Reeb graphs; 

• show that this distance function is stable under perturbations of the input data; 

• define ‘smoothing’ operations on Reeb graphs (which reduce topological complexity). 

The novelty of this paper is that we address these geometric questions using methods from 
category theory. Reeb graphs can be identified with a particular kind of cosheaf [24,48], and 
these cosheaves may be compared using an interleaving distance of the kind studied in [9,11]. We 
pull back that interleaving distance to obtain a distance function between Reeb graphs. While 
an efficient algorithm for computation of the interleaving distance is not yet available, one step 
of the process for construction yields a smoothed version of the given Reeb graph. This arises 
from a natural operation on cosheaves but has the added feature that it can be interpreted 
geometrically. We provide an explicit algorithm for constructing the smoothing of a given Reeb 
graph. The question of simplifying a Reeb graph, perhaps to deal with noise, has arisen in several 
applications [4,20,25,34], Our work differs from the solutions in those papers in that rather than 
doing local operations to collapse small loops, the smoothing operation is conducted globally and 
causes small modifications everywhere with the outcome that small loops are removed. 

Recently, other approaches to defining a metric between Reeb graphs have been defined; 
one is based on the Gromov-Hausdorff distance [4] and the other is defined using combinatorial 
edits [19]. These methods perhaps appear more natural from the geometric perspective, whereas 
our method is more natural from the sheaf-theoretic perspective. Our ideas have been inspired, 
in part, by the use of interleaving distances to compare join- and split-trees [30]. Finally, we 
point out that some of the category theory (without the cosheaves) appears in [41]; and a very 
extensive and accessible study of cosheaves can be found in [17]. 

Antecedents. It is well known amongst sheaf theorists that a locally constant set-valued 
cosheaf over a manifold M is equivalent to a covering space of M via its display space ; see Funk [24], 
Robert MacPherson observed that if a set-valued cosheaf on M is constructible with respect to 
a stratification (i.e. if it is locally constant on each stratum), then the cosheaf is equivalent to 
a stratified covering of M; see Treumann [43], Woolf [48] and Curry [17] for details. A stratified 
covering of the real line is what we call a Reeb graph. This leads to an equivalence between the 
category of Reeb graphs and the category of constructible cosheaves over the real line. 

Our definition of the interleaving distance between Reeb graphs is based on a very general 
framework for topological persistence developed by Bubenik and Scott [9] that was in turn 
inspired by the work of Chazal et al. [11] on algebraic persistence modules. Cosheaves are a 
particular kind of functor, as are generalized persistence modules, and the two ways of thinking 
overlap sufficiently to give us a metric on the category of Reeb graphs. 

In this paper, we give a quite detailed exposition of the ideas involved. We describe the 
equivalence of categories of Funk [24] explicitly in the situation that we need it, since the eventual 
goal is to use this equivalence in computations. Finally, whereas most of our work can be thought 
of as a combination of existing ideas from two separate fields, the smoothing operators we define 
on Reeb graphs seem to be novel, and hint at a richer family of operations on cosheaves to be 
discovered. 
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Figure 1: The Reeb graph is used to study connected components of levelsets. 

1.2 Reeb graphs and Reeb cosheaves 

Our starting point is a topological space X equipped with a continuous real-valued function 
/ : X —>■ [R. We call the pair (X, /) a ‘space fibered over R’ or, more succinctly, an R-space. 
For reasons of convenience we will often abbreviate (X, /) simply to /. The context will indicate 
whether we are thinking of / as a function or as an R-space. 

We can think of an R-space as a 1-parameter family of topological spaces f~ 1 (a), the levelsets 
of /. The topology on X gives information on how these spaces relate to each other. For instance, 
each levelset can be partitioned into connected components. How can we track these components 
as the parameter a varies? An answer is provided by the Reeb graph. 

The (geometric) Reeb graph of an R-space / is an R-space / defined as follows. First, we 
define an equivalence relation on the domain of / by saying two points i,i'eX are equivalent if 
they he on the same levelset / _1 (a) and on the same component of that levelset. Let Xj be the 
quotient space defined by this equivalence relation, and let / : Xj —y R be the function inherited 
from /. This is the Reeb graph. See, for example, Figure 1. 

If / is a Morse function on a compact manifold, or a piecewise linear function on a compact 
polyhedron, then its Reeb graph is topologically a finite graph with vertices at each critical value 
of /. This situation is well studied. These examples are included in a larger class, the constructible 
R-spaces, which have similar good behavior. We will say more about this in Section 2. If we work 
in greater generality, the quotient X —y Xj can be badly behaved. Among other things, we would 
need to pay attention to the distinction between connected components and path components. 
This is not an issue for constructible R-spaces, where the two concepts lead to the same outcome. 

We now indicate an alternate way of recording the information stored in the geometric Reeb 
graph. The abstract Reeb graph or Reeb cosheaf of an R-space / is defined to be the 
following collection of data (see Figure 2): 

• for each open interval ICR, let F(J) be the set of path-components of / _1 (J); 

• for / C J, let F [I C J] be the map F(J) —> F(J) induced by the inclusion / _1 (J) C / -1 (J). 

Let F denote the entirety of this data. It is easily confirmed that F is a functor (see Section 1.3) 
from the category of open intervals to the category of sets. As such, F is sometimes called a 
pre-cosheaf on the real line in the category of sets. The important point is that this information, 
in the constructible case, is enough to recover the geometric Reeb graph; see Figure 3. The 
other important point is that it is sometimes easier to work with the pre-cosheaf than with the 
geometric Reeb graph. 
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Figure 2: The geometric structure of a Reeb graph (left) may be represented categorically as a 
Reeb cosheaf (right). To each open interval we associate the set of connected components over 
that interval; to each inclusion of intervals we associate the map defined by component inclusions. 


There is considerable redundancy in the information stored by the sets F(J) and functions 
F[J C J] in the abstract Reeb graph of a constructiblc R-space. For example, the components 
over an interval / U J can be determined by considering the components over / and J and how 
they are related through the components over I D J. There are similar redundancies for every 
cover of an interval by other intervals. When systematized, these redundancies take on a standard 
form: they are precisely the conditions that ensure that the pre-cosheaf is a cosheaf. Thus, the 
abstract Reeb graph is renamed the Reeb cosheaf. 

We explain these standard ideas from sheaf theory more formally in Section 3. First we recall 
a few concepts from category theory, which provides the language for discussing these matters. 

1.3 Category theory 

We summarize what we need from category theory. For a general reference, see [29]. 

A category A is a collection of objects A E A, a collection of morphisms or arrows f : A A' 
between objects, and a composition operator that takes any two morphisms / : A —> A' and 
g : A 1 —> A" to a third morphism gf : A —> A" . The composition operator is associative and there 
is an identity morphism 1a ■ A —y A at each object A. There are many examples of categories 
found in all branches of mathematics. Here are some common examples: 

Category Objects Morphisms 

Set Sets Functions 

Vect Vector spaces Linear maps 

Top Topological spaces Continuous maps 

These are large categories, where the collection of objects is not a set but a proper class. 

Example 1.1: Any partially ordered set (P, <) can be thought of as a category P. The objects 
are the elements of P, and there is one morphism p —> q whenever p < q and no morphism 
otherwise. This is a small category, where the collection of objects is a set. 
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Figure 3: The Reeb cosheaf carries the same information as the Reeb graph. Here is one way 
the graph may be retrieved from the cosheaf. In the example, six carefully chosen intervals 
and their five pairwise intersections give rise to eleven sets (of components) and ten maps (of 
component inclusions). The graph built from this data—with eleven sets of vertices and ten sets 
of edges—has the same structure as the original Reeb graph. 



A functor F : A —> B is a map between two categories. It takes each object A 6 A to an 
object F(A) G B, and each morphism / : A —> A 1 of A to a morphism F[/] : F(A) —> F(A') of B, 
preserving composition and identities. A special case is the identity functor 1a : A —>■ A which 
takes each object and morphism to itself. 

A natural transformation ij : F G is a map between two functors F,G:A4B. It consists 
of a collection of morphisms r]A '■ F(A) —s- G(A), one for each object A 6 A, such that for each 
morphism / : A —> A' in A, the following diagram commutes: 


F(A) F(A') 


VA 


G (A) 


G[/] 


Va' 

G(A'). 


( 1 . 2 ) 


For any functor F, there is an identity natural transformation Ip : F F, defined at each 
object by (If) a = If (A)- Natural transformations can be composed in many different ways. In 
particular, if t] : F => G and 9 : G =>- FI then there is a composite Or] : F =>- FI defined at each 
object by (6*7/)^ = 9a^]a- These observations lead to the next example. 

Example 1.3: Let A, B be categories and suppose that A is small. Then the functors A —>• B 
themselves form a category B A , with natural transformations as the morphisms. 

A natural transformation 7] : F G is a natural isomorphism if each 7Ja is an isomorphism. 
By defining ( t]~ 1 )a = (va ) 1 we obtain the inverse natural transformation which satisfies 
Tj ^r] = If and = 1 g- Thus, natural isomorphisms are precisely the invertible morphisms 
in the functor category. 

Remark 1.4 (font convention): Some of the categories in this paper—specifically, Pre, Csh 
and Csh c —are categories of functors. The objects of these categories will be written in sans-serif 
style: F, G. We think of these as ‘small’ functors, the font style reminding us that we sometimes 
regard them as objects in a functor category. We contrast these with various ‘large’ functors 
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Figure 4: Road map of categories and functors. Categories of geometric objects (Section 2) 
occupy the left-hand column; categories of functors (Section 3) occupy the right-hand column. 
The up-arrows are inclusions of categories. The bottom row is an equivalence of categories. 


that are defined between the major categories of interest. These we write in calligraphic style: 

T,Q. 

It is often convenient to have more than one equivalent categorical interpretation of a given 
idea. 

• Two functors J-, Q are ‘essentially the same’ if there is a natural isomorphism between 
them, and we write J- ~ Q. Very often the isomorphism is canonically specified. This is 
much more common than the two functors being exactly equal to each other, which would 
be written T = Q. 

• Two categories A, B are ‘essentially the same’ if they are equivalent. This means that 
there is a pair of functors T : A — y B and Q : B —y A and a pair of natural isomorphisms 
H : QT =>■ 1a and v : TQ => 1b- An equivalence of categories refers to either the complete 
data ( } c , Q, fi,u) or one of the functors IF, Q by itself. 

Here we are mostly thinking of ‘large’ functors, as the font style suggests. 

1.4 Road map of categories and functors 

To develop the relationship between geometric and abstract Reeb graphs, we make use of several 
categories and functors. The reader may find it helpful to consult the road map in Figure 4. 
We define the various categories and functors in Sections 2 and 3, and establish the following 
relations: 

• 1ZT is naturally isomorphic to the identity functor on Reeb (Proposition 2.13). 

• C"1Z is naturally isomorphic to C (Theorem 3.17). 

• The functors C",T> define an equivalence between Reeb and Csh c (Theorem 3.22). 

Subsequently, we will define a metric on Pre and smoothing operators on R-Top and Pre. 
Through the diagram, these lead to a metric and a smoothing operator on Reeb. 
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2 The geometric categories 

In Sections 2.1, 2.2 and 2.3 we describe the three geometric categories, from largest to smallest. 
In Section 2.4 we define and study the geometric Reeb functor 1Z. 

2.1 The category of [R-spaces 

An object of R-Top is a topological space X equipped with a continuous map / : X —* R, denoted 
(X,/) or simply /. Point-preimages / _1 (a) are known as levelsets or fibers of the R-space. A 
morphism <p : (X, /) —>- (¥, g) is a continuous map ip : X —y Y such that the following diagram 
commutes: 



Composition and identity maps are defined in the obvious way. 

Remark 2.1: Being an example of a slice category , R-Top is sometimes named (Top j, R). 
In [41] it is called the category of scalar fields. 

2.2 The category of constructible [R-spaces 

We restrict to this class of spaces because the geometric Reeb graph of a general R-space may be 
badly behaved. These spaces are compact and have finitely many ‘critical points’ between which 
they have cylindrical structure. 

Formally, an object of R-Top c is an R-space that is isomorphic to some (X, /) constructed in 
the following way. A finite set of ‘critical values’ S = {cio, cq,..., a n } (listed in increasing order) 
is given. Then: 

• For 0 < i < n, we specify a locally path-connected compact space Vj. 

• For 0<i<n—1, we specify a locally path-connected compact space Ej. 

• For 0<z<n—1, we specify continuous maps lj : Ej —* Vj and iy : Ej —>• V i+1 . 

Let X be the quotient space obtained from the disjoint union of the spaces Vj x {cy} and Ej x 
[aj,aj + i] by making the identifications (lj(x),aj) ~ (ay ey) and (rj(x), a i+ i) ~ (ayay+i) for all i 
and all x G Ej. Let / : X —> R be the projection onto the second factor. 

Morphisms in R-Top c are the same as morphisms in R-Top (it is a full subcategory ). 

Example 2.2: The following R-spaces belong to R-Top c : (i) X is a compact differentiable 
manifold and / is a Morse function; (ii) X is a compact polyhedron and / is a piecewise-linear 
map; (iii) X is a compact semialgebraic subset of R n x R and / is the projection onto the 
second factor, (iii') X is a compact subset of R ra x R definable with respect to some o-minimal 
structure [16,44] and / is the projection onto the second factor. 

See Figure 5 for a manifold with a Morse function presented as a constructible R-space. 

Remark 2.3: The critical set is not uniquely specified, since we can always add extra critical 
points by splitting the cylinders Ej x [ay, a i+ i\ appropriately. One can request a minimal critical 
set, but we never specifically need it. 

The content of the next lemma is geometrically straightforward. We state it formally because 
we use it repeatedly to establish relationships between Reeb graphs and Reeb cosheaves. 




E 0 x [a 0 , ai] Eix[ai,a 2 ] E 2 x[a 2 ,a 3 ] E 3 x[a 3 ,a 4 ] E 4 x[a 4 ,a 5 ] 



¥ 0 Vi V 2 ¥ 3 ¥ 4 ¥ 5 


Figure 5: An example of a constructible R-space. Each of the six critical values of the Morse 
function on this torus contributes a critical fiber Vj, while the intervals between critical values 
contribute cylinders Ej x [a*, a i+ 1 ] which are attached to the critical fibers using maps defined by 
gradient flow. We encourage the reader to visualize the attaching maps explicitly. 


Lemma 2.4 (cylinder principle): Let (X, /) be constructible with critical set S = {ao, a 4 ,..., a n }. 
The fiber-inclusion maps 


Vj—» / 1 (a i - 1 ,a i+ i); x^(x,ai) 

Ej —» f~ 1 (a i ,a i+ 1 ); x 1 —y (x,a) some a G (a*, a i+ i) 

are homotopy equivalences which fit into diagrams 


= “ “ (2.5) 

Y ^ w 

f (fli— 1 ; ®i) ^ f (®j— 1 , ®i+i) f [flit ®i+i) 


that commute up to homotopy (0 < i < n, with E 0 , E n interpreted as empty spaces). The 
homotopy equivalences are natural, in the sense that we have commutative diagrams 


-»- / 1 (a i _i,a i+ i) 


E{ ->■ f 1 (a i ,a i+ 1 ) 


a 


a 

Y 


and 


a 


« ( 2 . 6 ) 

Y 


V 9 

l 


9 (®i—1; ®i+l) 


Ef 


9 (®i; ®i+l) 


whenever a : (X,/) —y (Y, g) is a morphism between R-spaces with critical set S. (For the left- 
hand maps we are identifying the spaces Vj, Ej with the corresponding fibers at eij, a respectively.) 

Proof. Thanks to the cylindrical structure between critical points, / _1 (aj,aj + i) is homotopy 
equivalent to any of its fibers Ej, and / _1 (aj_ 1 , aj +4 ) deformation-retracts onto its critical fiber Vj. 
The remaining assertions follow easily. □ 
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Eq x [ao, ai] £’ix[ai,a 2 ] E 2 x[a 2 ,a 3 ] E 3 x[a 3 ,a 4 \ £ 4 x[a 4 ,a 5 ] 



Vo Vi V 2 V 3 V 4 V 5 

Figure 6: An example of an IR-graph with its presentation as a constructive R-space. It is the 
Reeb graph of the R-space in Figure 5, its presentation obtained by applying the 7r 0 functor to 
the presentation indicated there. Indeed, the reader may verify that V* = 7r 0 (Vj) and E t = 7T 0 (Ej) 
and that the attaching maps are obtained in the manner described. 

2.3 The category of [R-graphs 

An object of Reeb, also known as an R-graph, is a constructible R-space (X, /) for which the 
spaces Vj and E, are O-dimensional (i.e. finite sets of points with the discrete topology). Geo¬ 
metrically, it is a compact 1-dimensional polyhedron triangulated so that restriction / to each 
edge E C X is an embedding. Morphisms in Reeb are the same as morphisms in R-Top. 

Notation 2.7: We construct an R-graph (X, /) with critical set S = {ao, ai,..., a n } as follows: 

• For 0 < i < n, we specify a finite set of vertices V,, which lie over a*. 

• For 0 < i < n — 1, we specify a finite set of edges Ei which lie over the interval [a*, Oj+i]. 

• For 0 < i < n — 1, we specify attaching maps £i : Ei —> V) and r* : Ei —» V l+ \. 

The space X is the quotient of the disjoint union of the spaces Vj x {a*} and E t x [a*, a i+ 1 ] with 
respect to the identihcations (^( e ); a «) ~ (e,aj) and (rj(e), aj+i) ~ (e,a, + i), with the map / 
being the projection onto the second factor. See Figure 6. If we wish to add extra points to the 
critical set, we can retain this description by splitting the edges at new vertices over the new 
critical points. 

The restriction of / over each open interval (aj,aj + i) is a covering map. We use this fact in 
the next proposition, which gives a combinatorial description of the morphisms of Reeb. 

Proposition 2.8: Let (X, /), (¥, g ) be R-graphs with a common critical set S = {ao, a 4 ,..., a n } 
and described as above. A morphism <p : (X,/) —* (Y, g) is exactly specified by the following 
data: 

• Maps iff ■ V/ —* Vf for 0 < i < n. 
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• Maps ipf : E- —>• Ef for 0 < i < n — 1. 

• The consistency conditions gffd{ = tf<pf and <pY+i r I = uVf are satisfied, for 0 < i < n — 1. 

Proof. Any morphism <p : (X, /) —> (Y, g) defines consistent vertex and edge maps as above (the 
covering map structure between critical points guarantees that each edge of X maps to exactly 
one edge of ¥, in a unique way once the edge is chosen). Conversely, a morphism can be specified 
by defining continuous maps on the vertices and edges of X in a consistent way; the data above 
provide that. □ 

The requirement of a common critical set is no restriction, because we can take the union 
of critical sets for (X, /) and (Y, g) to obtain a critical set for both. On the other hand, for 
computational purposes we may wish to be more economical; see Section 5.3. 


2.4 The Reeb functor 7 Z 


The Reeb functor 1Z converts a constructible R-space to an R-graph, its Reeb graph. We 
provisionally define it as a functor R-Top —>■ R-Top, then show that it restricts to a functor 

R-Top c —>Reeb. 

Lemma 2.9 (quotient principle [40, Proposition 3.8.2]): Let (X/~) be the quotient of a topo¬ 
logical space X by an equivalence relation and let Y be another topological space. For any con¬ 
tinuous function X —y Y which is constant on equivalence classes, the induced function (X/~) —y Y 
is continuous with respect to the quotient topology. □ 

Let (X,/) be an R-space, with geometric Reeb graph (Xy,/). Recall that Xy is the quotient 
space of X by the relation whose equivalence classes are the path-components of the levelsets 
of /. It follows from the quotient principle that / is continuous, so (X/,/) is an R-space. The 
quotient map X —)• X/ defines a morphism (X, /) —> (X/, /) in the category of R-spaces; we label 
this morphism Pf = P(x,f)- 

Now consider a morphism a : (X,/) — > (Y, g). Since a preserves levelsets and (being con¬ 
tinuous) carries path-connected sets to path-connected sets, the composite map X —)• Y —> Y g 
is constant on equivalence classes. By the quotient principle, the induced map a : X/ —* Y g is 
continuous. Since ga = f this defines a morphism a : (X/, /) —>■ (Y g ,g). 

Proposition 2.10: The formulas 72.(X, /) = (X/, f) and lZ[a] = a define a functor R-Top —>■ 
R-Top. The collection of maps p = (p/) constitute a natural transformation 1r_t op => 1Z. 

Observation 2.11: In other words 1Z is a pointed endofunctor of R-Top: an endofunctor that 
is the target of a natural transformation from the identity functor. We call 1Z the Reeb functor 
and its ‘basepoint’ p the canonical projection. 


Proof of Proposition 2.10. Notice that a is the unique map that makes the following diagram 
commute: 


X —^-3- y 


Pf 

v 


X/ 


Pg 


a. 

- 3 ^ 




Uniqueness implies that 1Z[] respects identities and composition, so 1Z is a functor. Of course, 
these facts are easily verified directly. The commuting of the square is what makes p a natural 
transformation. D 


Proposition 2.12: The Reeb functor carries constructible R-spaces to R-graphs. 
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To see how this works, compare Figure 5 to Figure 6. 

Proof. Given a constructible R-space (X, /) it is clear how 7£(X, /) should be described as an 
R-graph (using Notation 2.7): 


V = 7To(Vi), = 7r 0 (Ei), ^ = 7T 0 [li], r* = 7r 0 [r ? ;]. 

We have a tautological bijection from X/ to this graph G, since the points of G are precisely the 
path-components of the levelsets of /. This map preserves levelsets. It remains to show that it 
is a homeomorphism. 

First we show that it is continuous. Let X denote the disjoint union of the spaces Vj x {a,} 
and E,; x [a,i,a i+ 1 ]. Consider the composite X ^ X Xj -)■ G; the first two maps are quotient 
maps. This composite is continuous because the path-components of a locally path-connected 
space are open. Applying the quotient principle twice, it follows that X/ —* G is continuous. 

The proof is completed by invoking the standard result that a continuous bijection from a 
compact space (Xj) to a Hausdorff space (G) is a homeomorphism [40, Theorem 5.9.1]. □ 

Proposition 2.13: Each R-graph is naturally isomorphic to its Reeb graph. 

Proof. Each levelset of an R-graph is a finite discrete space, so the equivalence classes are sin¬ 
gletons. Thus the canonical projection p is a homeomorphism in these cases. O 

Thus we can think of the Reeb functor as a projection operator R-Top c —y Reeb. Henceforth, 
we will mostly reserve the symbol 1Z for the functor with this particular domain and codomain. 
Proposition 2.13 can be restated as the assertion that p restricts to a natural isomorphism 

lReeb => 77Z. 

3 The cosheaf categories 

We now describe the three cosheaf categories, from largest to smallest, in Sections 3.1, 3.2 and 3.3. 
In Section 3.4 we study the Reeb cosheaf functor C, and in Section 3.5 we show that it defines 
an equivalence of categories. 

The idea behind the cosheaf categories is that we can study an R-space by inspecting its 
behavior over subintervals of R. To this end, let Int denote the category whose objects are 
open intervals ICR and whose morphisms are the inclusions / C J. (This is an instance of 
Example 1.1.) 

3.1 The category of pre-cosheaves 

The largest of the three cosheaf categories is Pre = Set Int , the category of functors Int —» Set 
with natural transformations as morphisms (Example 1.3). The elements of Pre are called 
pre-cosheaves. 

Remark 3.1: More generally for any category C we can consider Pre(C) = C Int , the category 
of pre-cosheaves in C over the real line. 

Example 3.2: Let (X, /) be an R-space. This determines a pre-cosheaf X e Pre(Top) over 
the real line, as follows: for every interval / we let X(J) be the topological space / _1 (J), and for 
every pair I C J we let X[J C J] be the inclusion map / _1 (J) C / _1 (J). 

Example 3.3: The previous example generates many others. We can post-compose X with 
any functor Q : Top —y C to obtain a pre-cosheaf QX G Pre(C). For example: 
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• Let Hfc denote singular fc-homology; then H^X is a pre-cosheaf in Ab, abelian groups. 

• Let 7 Tq denote the set of path-components of a space; then 7r 0 X is a pre-cosheaf in Set. 

• Let 7fo denote the set of connected components of a space; then tTqX is a pre-cosheaf in Set. 

Thus 7 r 0 X(J) = 7r 0 (/ _1 (/)), for example. The key requirement is that the operations H fc , 7r 0 , 7f 0 be 
functors: they specify an object for each topological space, and a morphism for each continuous 
map. For instance, if p : Y —>■ Z is a continuous map then each path-component of Y maps into 
a path-component of Z. This defines 7To(y] : 7To(Y) —> 7T 0 (Z). 

3.2 The category of cosheaves 

The second category in the right-hand column is Csh, the category of cosheaves in Set over 
the real line. It is a full subcategory of Pre: it is defined by specifying which pre-cosheaves are 
cosheaves, and declaring that cosheaf morphisms are the same as pre-cosheaf morphisms. 

A cosheaf is a pre-cosheaf F which satisfies the following ‘gluing’ property. Let U be an open 
interval and let (I p \ p € P) be a family of open intervals whose union is U. Then we ask that 
F (U) be the colimit of the following diagram: 

H F(/„ n I,) =t n F(/„) (3.4) 

p,q v 

This must be true for every U and every cover (/ p ). In particular, this implies F(0) = 0 so on 
the left-hand side of (3.4) we consider only the terms with I p fl I q nonempty. 

Here are three interpretations of the gluing property. 

1. F (U) is obtained from the disjoint union ]j p F(/ p ) by identifying all pairs of points 

F [i p n i q c ip](x) e F (ip) 

F[/ p n4c/,](i)eF(/,) 

where p, q are indices with I p fl I q nonempty and x G F(/ p fl If). 

2. F (U) is the set of connected components of the graph with a vertex for every element in the 
disjoint union F(/ p ) and an edge for every element in the disjoint union q F(/ p fl If). 
The maps F [I p fl I q C I p ] and F[/ P fl I q C If indicate the vertices to which each edge is 
glued. 

3. F (U) is characterized by the following universal property. Let Z be a set, and suppose that 
maps ( p : F(J p ) —> Z are given for all p, such that for all p, q with I p fl I q nonempty the 
following two maps F(/ p fl If) Z are equal: 

C V ° n Iq ^ 4] = Cq ° F [I p fl I q C If (3.5) 

Then there is a unique map ( : F([7) —> Z such that 

C P = C°F[J p Cf/] (3.6) 


for all p. 

The first two interpretations are valid in Set. The third interpretation is meaningful in any 
category. The universal property characterizes F (U), and the maps to it from the F(/ p ) and the 
F (4 fl I q ), uniquely up to a canonical isomorphism. 
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Example 3.7: In Example 3.3 above, II/X is not in general a cosheaf. Consider open intervals 
U = I U J and N — I D J. The Mayer-Vietoris theorem gives the following exact sequence of 
abelian groups: 

H*(/-'(JV)) —>■ H, .(/-'(/)) e A H t _ 1 (/- 1 (JV)) —►... 

If the map labeled d were zero, then the exactness of this sequence implies that the gluing 
condition holds for the cover U = / U J. When d is not zero, the gluing condition fails for this 
cover. 

Example 3.8: Continuing with Example 3.3, we will later see that 7ToX is always a cosheaf 
(Proposition 3.15), but 7f 0 X sometimes is not (Example 3.20). 


3.3 The category of constructible cosheaves 

The third category in the right-hand column is Csh c , the category of constructible cosheaves 
in Set over the real line. It is a full subcategory of the category of cosheaves, defined by specifying 
which cosheaves are constructible and using the same morphisms as before. 

Definition 3.9: A cosheaf or pre-cosheaf F is constructible if each F(J) is finite and there exists 
a finite set S' C [R of ‘critical values’ such that: 

• if / C J are open intervals with / fl S — J D S then F[J C J] is an isomorphism; 

• if / is contained in (—oo, min(S')) or (rnax(S'), +oo) then F(J) is empty. 

As with constructible R-spaces, if the conditions hold for some S then they hold for any S' D S. 

For a given critical set S = {cio, Gq, ..., a n } the following ‘zigzag’ diagram in Int is of particular 
importance: 


(—oo, a\) 



(ao,cq) 




(do, Oi) 


(ai, (X2) 


ifln—2i ®n) 


(dn—1, +00) 




(®n— 1 


O'n) 


(3.10) 


Notation 3.11: For a constructible cosheaf F with critical set S, write 

V[ = F((dj_i, d i+ i)), for 0 < i < n 
= F((dj, d i+ i)) for 0 < i < n — 1 

where d_i = —00 and a n+ i = +00; and 

A F = F[( a i) a i+l) ( a i-l) a i+l)] : El —> V) F 
rf = F[(di, d i+ i) C (di, a i+2 )\ : Ef V^ +1 


for 0 < i < n — 1. 

As we will see in Proposition 3.18, this collection of combinatorial data suffices to define a 
constructible cosheaf that is unique up to canonical isomorphism. To get to this result—and 
more importantly to move towards the equivalence of categories, Theorem 3.22— we establish a 
combinatorial description of morphisms between constructible cosheaves. 
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Proposition 3.12 (combinatorial description of cosheaf morphisms): Let F, G be constructible 
cosheaves, and let S = {ao, cq,..., a n } be a common critical set for them. A morphism 0 : F =>• G 
gives rise to the following data: 

• Maps 00 : Vf —» I0 G for 0 < i < n. 

• Maps 00 : Ef —* Ef for 0 < i < n — 1. 

• Consistency conditions: 00 0 F = 7 G 00 and 00_ 1 rf = r f00 for 0 < i < n — 1. 

Conversely, any collection of maps 00,00 satisfying the consistency conditions arises in this way 
from a unique morphism 0 : F =>• G. 

Remark 3.13: As usual, the requirement of a common critical set is no restriction, because 
we can take the union of critical sets for F and G to obtain a critical set for both. 

Proof. The first assertion is clear: we set 00 = ' l l ) {a i - 1 ,a i+1 ) an d 'iff = and the consistency 

conditions follow from the naturality of if with respect to the inclusions {a i: a i+ 1 ) C (aj_i,a i+ i) 
and (cij, ^ (®i, cq+ 2 )• 

For the converse, we must show that the collection of maps 0( ai _i,a i+ i) = 0f and 'if(< H ,a i+1 ) = 
iff extends uniquely to a natural transformation 0 = (-0/ | / E Int). 

We begin by showing that 0/ is uniquely determined for intervals of the form I = (Opa*,). 
The starting data provides these maps for intervals of ‘length’ (i.e. k — j equal to) 1 or 2. Longer 
intervals can be expressed as a union of intervals of length two. We then have: 

' k—2 k—2 

F((ai,a i+ i)) JjF((ai,a i+2 )) 

.i=j+l i=j 

' k—2 k—2 

J J G((a*, a i+ 1 )) =1 ]J G((ap a i+2 )) 

i=j+1 i=j 

By naturality, the desired map 0 ( a . ak ) must be compatible with the maps 0^,0® that map the 
terms in the colimit for F to the terms in the colimit for G. I 11 particular, 0( aji a fe ) must factor the 
maps 

Ci = G[(aj, a i+2 ) C (aj,a k )] o 0f : F((ai,a i+2 )) ->■ G((a i ,a fc )). 

The consistency conditions imply that the family (0) satisfies (3.5). From the universal property 
for colimits, there is a unique 0( a ,a fc ) = C compatible with the 00,00. 

For an inclusion I C J between two such intervals, the naturality condition is that the square 


F(/) - 

-^.F (J) 



G (/) ■ 

GKJ| - G(J) 


F a k )) = colim 

G ((dj, ah)) = colim 


commutes. If / is an interval of length 1 or 2, this is precisely the compatibility demanded in the 
construction of 0j. If / is a longer interval, then the diagram commutes when F(J) is replaced by 
any of the the individual terms in its colimit diagram. The universal property for colimits then 
implies that the two sides of the square F(J) —> G( J) are equal, being the unique map compatible 
with the map on each individual term. 

To finish, we consider arbitrary open intervals. Any such I is contained in a unique maximal 
interval / = ( aj,a k ) that meets S' in the same subset. Since G[J C I] is an isomorphism, we can 
and must define 

0/ = G[J C 7]- 1 o0 | oF[7 C I] 
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to satisfy naturality for / C I. That done, naturality for / C J follows from naturality for 
J C J. □ 

Corollary 3.14: Let ip : F => G be a morphism between constructible cosheaves with common 
critical set S = {ao, cq,..., a n }. If xpi is an isomorphism for each ‘short interval’ / = (cq_i, di + 1 ) 
and / = (ai,di + i) then ^ is a natural isomorphism; that is, an isomorphism of cosheaves. 

Proof. We construct the inverse u : G => F by setting uii = ipj 1 for short intervals / and 
applying the proposition (consistency follows from the naturality of ip). Since (c uip)i = 1f( 7) and 
{ipuS)i = 1g(j) for short intervals, the proposition implies ujip = Ip and ipu = 1g- Thus ip is a 
natural isomorphism. □ 


3.4 The Reeb cosheaf functor C 


The Reeb cosheaf functor, C, converts an R-space to its Reeb cosheaf. Let / = (X, /) be an 
R-space. Then C(f) — F is the pre-cosheaf defined by 

F(7) = n „/-'(/), F[7 C J] = Jr 0 [/ _1 (/) C / -1 (J)]. 

(This is the second item in Example 3.3.) 

Any morphism a : (X,/) — y (Y ,g) yields a natural transformation of pre-cosheaves C[a] : 
F =>■ G defined as follows. For each interval /, the map a carries /~ 1 (J) into g -1 (J). We set 

C[a]i = 7r 0 [/ _1 (/) A ( 7 ^ 1 (/)] which is a map F(J) — > G(/). 


We prove the naturality of C[a] by applying 7T 0 to the commutative square on the left: 


rv) 


9~V) 


r\j) 


9~\J) 


WO, 


F(/f-^U F(J) 


C[a]/ C[a], 

G (I) G(J) 


Since C[—\ respects composition and identities, we have a functor from R-spaces to pre-cosheaves. 

Proposition 3.15: The pre-cosheaf F = C(f) is a cosheaf. 

Proof. Let U be an open interval and let (I p \ p G A) be a cover of U by open intervals. We 
show that F (U) satisfies the universal property for the colimit of (3.4). Accordingly, let Z be a 
set and let ( p : F (I p ) —> Z be functions satisfying the consistency condition (3.5). We show that 
there is a unique ( : F(f7) —>• Z satisfying Eqn. (3.6). 

Let [y\u denote the path-component of a point y in / _1 (17). Since f(y) must belong to 
some I p , we are forced to define 

C(Mc/) = C ° F [Ip c U\{[y] Ip ) = Cp([y]i P ) 

and moreover the right-hand side does not depend on the choice of p, because 

Cp{[y]i P ) = C p° F [ip n/,c ip]{[y]i P rii q ) = C<? ° F [i p n i q c. iq]{[y]i p ni q ) = C q{[y]i q ) 
by condition (3.5), if f(y) £ I p fl I q . 

It remains to show that this definition is independent of the point y used to identify the 
component. Suppose [yf\u = [; yi]u ■ Then there is a continuous path (y t ) in f~ l {U) from y 0 to y\. 
Now every point in [0,1] has a neighbourhood over which f(y t ) is contained in some fixed I p . 
Over that neighbourhood, [yt]i p is constant and therefore Cp( [yt] r p ) is constant. Since [0,1] is 
connected, this local constancy implies global constancy and so C([2/o]t/) = C([2/i]c/)- □ 
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Proposition 3.16: If / is a constructible R-space then F = C(f) is a constructible cosheaf 
with the same critical set. 


Proof. Let / C J be intervals which meet S in the same set of points. The product structure over 
the components of R \ S' implies that / _1 (J) C / _1 (J) is a homotopy equivalence and therefore 
F[J C J] is an isomorphism. Q 

It follows from Propositions 3.15 and 3.16 that the operation C defines functors as follows: 

C : R-Top ->■ Csh, C : R-Top c ->■ Csh c , C" : Reeb ->■ Csh c . 

We use the symbols C' and C" when we wish to be precise about the domain and range of our 
functors. When that is not important, we simply write C. 

Theorem 3.17: The functors C' and C"1Z : R-Top c —> Csh c are naturally isomorphic. 

In other words when starting with a constructible R-space, we can immediately use C to 
convert it to a cosheaf or we can take its geometric Reeb graph and then use C to convert it to 
a cosheaf; either way the result is the same. 

Proof. First, we compare C, C1Z regarded as functors R-Top — y Csh. There is a natural transfor¬ 
mation rj = Cp defined by applying the functor C to the canonical projection p (Observation 2.11). 
Spechcally: 

ri, = C[p,} = C[(X, /) -> (X/, /)] : C(f) —> CR(/) = C(f) 

We must show that rj is an isomorphism at each constructible R-space, meaning that it restricts 
to a natural isomorphism C =>• C"1Z. Let / = (X,/) be constructible; then the cosheaves C(f), 
ClZ(f) are themselves constructible with the same critical set. To show that rjf : C(f) => ClZ(f ) 
is an isomorphism it is enough, by Corollary 3.14, to show that (rjf)i : [C(/)](/) — > [ClZ(f)\(I) is 
an isomorphism whenever / is a short interval. 

The cases I = (aj_i, a i+ i) are handled by the left diagram, and the cases / = (aj,a i+ 1 ) are 
handled by the right diagram: 

V| ^ r\l) E, f~\l) 


| Y Y 

v; r\i) Ei /-'(/) 


This is (2.6) from the cylinder principle (Lemma 2.4) for the canonical projection (X, /) —> 
(X/,/). The horizontal inclusions are homotopy equivalences, and the left-hand map of each 
diagram induces a bijection of path-components, so the same is true of the right-hand maps. In 
each case, we see that (rjf)i = 7r 0 [/~ 1 (/) —> / _1 (/)] is an isomorphism. □ 

We round out this section with a result promised earlier and two cautionary examples. 


Proposition 3.18 (combinatorial description of constructible cosheaves): Given a critical set 
S = {a 0 , ai,..., a n }, finite sets 


V 0 ,...,V n and E 0 ,...,E n _ i 

and maps 

7* : Ei —>■ Vi and r^ : Ei —> R+i 

for all i. Then there exists a constructible cosheaf F with critical set S, together with (in 
Notation 3.11) bijections Vf = Vj and Ef = Ei such that the maps and rf correspond to the 
maps £i and r t . Any two such cosheaves are canonically isomorphic. 
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This fact, together with Proposition 3.12, amount to a particular instance of a more general 
result of MacPherson that we discuss in Section 6. 

Proof. We may construct F as the Reeb cosheaf of the IR-graph with critical set S constructed 
from the same combinatorial data. The conditions on F are easily verified; and it is a cosheaf with 
critical set S, by Propositions 3.15 and 3.16. Given another such cosheaf G, we have canonical 
identifications Vf = Vj = V G and Ef = Ei = Ef through which £f,rf correspond to £ t , r t and 
then to £f , rf . These identifications dehne an isomorphism of cosheaves, by Corollary 3.14. Q 

Remark 3.19: For readers more familiar with category theory, the existence of the cosheaf F 
may be proved more directly (i.e. without manufacturing a topological space) as follows. The 
starting data specifies the value of the cosheaf on every open interval that meets at most one 
critical point. Every interval I is the union of its subintervals of this type. We define F(J) to be 
the colimit associated to this union, and then invoke general properties of colimits to show that 
F is a cosheaf with the desired values on short intervals. 

The good behavior of our functors on constructible [R-spaces does not extend to the non- 
constructible case. Here are two counterexamples based on the topologist’s sine curve 

S = {(x,sin(l/x)) | 0 < x < 1} U {(0,2/) | y G [-1,1]} 

which is a connected but not path-connected compact subset of the plane. 

Example 3.20: Proposition 3.15 fails if we replace the path-component functor 7To by the 
connected-component functor 7f 0 . Consider (S, y) where y is the projection onto the second 
coordinate. Now S itself is connected. On the other hand, if / is an interval that meets but 
does not contain [—1,1] then / _1 (J) consists of countably many connected components, one of 
which is the segment on the y-axis. Now cover the real line by intervals of that type. The 
associated colimit has at least two elements, since the segment on the y-axis is always separate 
from everything else. This breaks the cosheaf condition, since 7fo/ _1 (lR) = 7fo(S) is a singleton. 

Example 3.21: The natural isomorphism of Theorem 3.17 does not extend to arbitrary R- 
spaces. Consider (S, x) where x is the projection onto the first coordinate. We can identify its 
Reeb graph as follows: each levelset over [0,1] is path-connected, so the projection to the x-axis 
induces a continuous bijection, and therefore homeomorphism, from the Reeb graph S x to the 
interval [0,1]. Then 7ro(S x ) = 7r 0 ([0,1]) is a singleton, whereas 7r 0 (S) is not. We deduce that the 
cosheaves C(S,x) and C77(S,x) = C(S x , x) are not isomorphic: they return non-isomorphic sets 
when evaluated at the interval R (or indeed any interval containing 0). 

3.5 Equivalence of categories 

This is the theorem that relates Reeb graphs to Reeb cosheaves. 

Theorem 3.22: The functor C" : Reeb —> Csh c is an equivalence of categories. 

Proof. We will show that C" = C|R ee b is fully faithful and essentially surjective. This means: 

(i) For every / = (X, /) and g = (¥, g ) in Reeb, the map 

Hom(/,g) ^GiHom(C(/),C( 9 )) 

is a bijection of sets. 

(ii) For every F G Csh c there exists / = (X, /) in Reeb such that C(f) is isomorphic to F. 
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It is a theorem [29, §IV.4] that such a functor is an equivalence of categories. 

(i) We show that C" is fully faithful. Let / G Reeb and write F = C(f). With respect to a 
critical set for both, we can describe / and F by data 

V /, E[ , i{ , r{ (Notation 2.7) 

1/ F , , l\, rf (Notation 3.11) 

respectively. Applying 7r 0 to the homotopy equivalences from the cylinder principle (Lemma 2.4) 
and to diagram (2.5), we get isomorphisms 

V/ = MV/) =n 0 f-\a i - 1 ,a i+1 ) =V? 

E{ = ME() =tt 0 r\a h a i+1 ) = E\ 


which carry i{,r{ to £f,rf. 

Now let /,g G Reeb and write F = C(/), G = C(g). With respect to a common critical 
set, we can describe f,g and F, G by corresponding sets of data, with isomorphisms as above. 
From the characterizations of morphisms given in Propositions 2.8 and 3.12, there is an obvious 
bijection between Hom(/, g) and Horn (C(f),C(g)) 

f V y : V/ -► V? \ I : VI -> Vf \ 

\ V f : E{ —^ Ef J ^ \ : Ef Sf / 

defined using these isomorphisms. To show that this bijection is given by C[—\ we apply 7r 0 to 
the diagrams (2.6) from the cylinder principle. This completes the proof that C" is fully faithful. 


(ii) We show that C" is essentially surjective. Let F G Csh c with critical set S = {ao, cq,..., a n }. 
Let (X, /) be the IR-graph with critical set S defined by the following data (Notation 2.7) 

Vi F((dj_i, cij_|_i)) F[(aj, a i+ i) C (cij_i, (Zj_)_i)] (3.23) 

Ei = F((aj,a i+ i)) r t = F[(aj,a i+ i) C (ai,a i+2 )\ (3.24) 

and write F 7 = C(f). 

From the cylinder principle (Lemma 2.4) and diagram (2.5) we obtain isomorphisms 

F , ((ai-i, Oj+i)) — Vi = F((aj_i,a i+ i)) (3.25) 

F / ((a*, aj+i)) = Ei = F((ai,a i+ i)) (3.26) 

which are natural with respect to the inclusions (aq aj+i) C (a*- 1 , a,: + i) and (a*, a* + i) C (a*, 0 , 1 + 2 ) ■ 

Then Corollary 3.14 implies that the cosheaves F and F 7 = C(f) are isomorphic. □ 

Since C" is an equivalence of categories, it has an inverse functor T> : Csh e —>• Reeb called 


the display locale functor [24], There are various ways to define this functor; the result is 
unique up to a canonical natural isomorphism. We can define T> combinatorially as follows: given 
F G Csh c with critical set S, let V( F) be the R-graph defined by the data (3.23) and (3.24). A 
natural isomorphism CV lcsh c is defined by the identifications in (3.25) and (3.26). These 
identifications uniquely determine the result T>[a\ of applying P to a morphism a. 

Corollary 3.27: The functors T>C', 1Z : [R-Top c —y Reeb are naturally isomorphic. 

That is, the Reeb graph of a constructible R-space is equal to the display locale of its Reeb 
cosheaf. 
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Figure 7: Road map with functors for smoothing cosheaves (S £ ), for thickening [R-spaces (71), 
and for smoothing R-graphs (U e ). 


Proof. We have 77 ~ VC'IZ ~ VC by Theorems 3.22 and 3.17. O 

Remark 3.28: The display locale of a cosheaf can be defined more abstractly and generally [24] , 
yielding a functor V : Csh — y R-Top. The fiber of 77(F) at a G R is defined to be the limit of 
F(J) over intervals I containing a. This is called the co-stalk of the cosheaf at a. The disjoint 
union of these co-stalks is topologized as follows: for each interval / and x G F(J) there is a basic 
open set Uj, x defined to be the elements of the co-stalks at all a G / which project to x. 

4 The interleaving distance 

We are ready to define the distance between a pair of Reeb graphs. The abstract principle is quite 
simple (Section 4.1): we regard the Reeb graphs as constructible cosheaves, then we compare 
the cosheaves using an ‘interleaving distance’ [9,11]. Of course, we would like to interpret this 
as geometrically as possible. To do this, we consider two parallel operations: smoothing of pre¬ 
cosheaves (Section 4.2) and thickening of R-spaces (Section 4.3). The interleaving distance may 
be expressed in terms of smoothings; the resulting distance on Reeb graphs may be expressed in 
terms of thickenings. 

The smoothing functors (S e ) and the thickening functors (T e ) give compatible transforma¬ 
tions on the two sides of our road map: see Figure 7. Smoothing preserves the subcategories 
of cosheaves and constructible cosheaves. In a sense, that explains why it has geometric signif¬ 
icance and why the existence of the thickening functor is not a surprise. Thickening preserves 
the subcategory of constructible R-spaces. This allows us to define a semigroup of topological 
smoothing functors (U e ) on Reeb graphs (Section 4.4). 

4.1 Interleaving of pre-cosheaves 

Interleavings of persistence modules were used, implicitly, in the proof of the persistence stability 
theorem of Cohen-Steiner et al. [13]. Chazal et al. defined the concept explicitly for their algebraic 
stability theorem [11]. More recently Bubenik and Scott have given a general formulation in 
categorical language [9]; we follow their ideas closely. 

Interleavings are approximate isomorphisms. Let F, G : Int —>• Set be pre-cosheaves. Recall 
that an isomorphism between them consists of two families of maps 

<Pi : F(/) G(J), Vh : G(/) -G F(J) 
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that are natural with respect to inclusions / C J, such that <p/,0/ are inverses for all /. 

We can give ourselves e leeway by expanding the intervals slightly. For any interval I = (a, b ), 
let I £ = (a — e, b + e) denote the interval expanded by £ > 0 . 

Definition 4.1: An £-interleaving between F, G is given by two families of maps 

<Pi : F(J) -> G(/ £ ), Vh : G(J) F(/ £ ) (4.2) 

which are natural with respect to inclusions / C J and such that 

o ip, = F[J C I 2e }, ^o^; = G[/C I 2e ] (4.3) 

for all I. When e = 0 this is is precisely an isomorphism between F, G. 

Definition 4.4: The interleaving distance between two co-presheaves F, G : Int —> Set is 

defined 

d;(F, G) = inf (e | there exists an £-interleaving between F, G). 

The Reeb distance between two IR-graphs / = (X, /) and g = ( Y,g ) is defined 

dn(f,g) = d l (C(f),C(g)). 

(The infimum of an empty set is understood to be oo.) 

This definition of dR may not seem immediately helpful: it requires converting the Reeb 
graphs into cosheaves, and then comparing the cosheaves by a metric that is itself somewhat 
mysterious. In the next few sections we will develop a more geometric formulation that allows 
us—at least in principle—to compute the distance function. Having said that, certain geometric 
assertions are immediately accessible. We present some of these results now. 

Proposition 4.5: The interleaving distance di is an extended pseudometric on Pre: it takes 
values in [0, oo], it is symmetric, it satisfies the triangle inequality, and di(F, F) = 0. It follows 
that dR is an extended pseudometric on Reeb. 

Proof. For the triangle inequality, note that if (^}), (fi}) define an £i-interleaving between F, G 
and {v>i)i (^/) define an £ 2 -interleaving between G, FI then 

v] = v]ci °v)i -0f = V / e 2 ° 0? 

define an (ey + £ 2 )-interleaving between F. H. The remaining statements are obvious. □ 

This approach to interleaving distances [9,11] is designed to make the next theorem as easy 
as possible. 

Theorem 4.6 (Stability of Reeb distance): (i) Let (X, /) and (X, g) be [R-spaces (with the same 
total space X). Then: 

di(C(/),C(g)) < II/ —fliloo 

(ii) Let (X,/) and (X, < 7 ) be constructible [R-spaces. Then the Reeb graphs (X/,/) and (X 9 ,g) 
satisfy: 

d r(/,<?) < ||/-Slice 

Proof, (i) Suppose ||/ — g||oo < e. We show that there is an £-interleaving between C(f),C(g). 
The supremum bound implies that we have inclusions 

f-\l)Qg~\l E ), g-\l) C f~\F) 
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for all /. Accordingly, we define 

<Pi = Mr\l) Q 9~\ ie )l = ko [g _1 (-0 Q r l {i £ )\- 

Naturality and the other conditions are satished because diagrams of inclusions always commute. 

(ii) This follows from part (i) because the natural isomorphism C"1Z = C (Theorem 3.17) 
implies that C(f),C(g ) are isomorphic to C(f),C(g). f3 

The Reeb distance is sometimes infinite. 

Proposition 4.7: The Reeb distance between two R-graphs (X,/), (Y ,g) is finite if and only 
if X, Y have the same number of path components. 

Proof. Let I be a very large interval containing /(X) U g(Y). Writing F = C(f) and G = C(g) we 
have 

7T 0 (X) = F(J) = F(/ e ) = F(/ 2e ), vro(Y) = G(J) = G(/ £ ) = G(/ 2e ) 

for any e > 0. Thus any e-interleaving defines a bijection 7r 0 (X) = 7r 0 (P) through its maps tpi, ifi. 

Conversely, suppose there is a bijection 7r 0 (X) = 7r 0 (Y). Let e be larger than the diameter of 
/(X) U g(Y). This implies that if / meets either /(X) or g(Y) then F(/ £ ) = 7T 0 (X) and G(/ e ) = 
7 To(Y). For these intervals we define (pi,ifi as the following composites, using the bijection. 

F[/Ci»] 

<Pi ■ F(/) -> F(/ e ) = 7T 0 (X) = 7T 0 (Y) = G(P) 

G[ICI S ] 

if! : G(/)-> G(/ e ) = vr 0 (Y) = 7T 0 (X) = F(/ £ ) 

For the remaining intervals, F(J) = G(J) = 0 so nothing needs to be done. It is not difficult to 
verify that these maps define an e-interleaving. □ 

Proposition 4.8: The Reeb distance between two R-graphs is zero if and only if they are 
isomorphic. 

Proof. The nontrivial part is to show that Reeb distance zero implies that the R-graphs are 
isomorphic. We will show that if F, G are constructible cosheaves with d;(F, G) = 0 then they are 
isomorphic. This is an equivalent statement since C" is an equivalence of categories. 

Here is a quantified statement that implies the result. Let S = {a 0 , oq,..., a n } be a common 

critical set for F, G and let H = min ? ;(aj + i — af) > 0. We claim that if F, G are e-interleaved for 

£ < h /4 then F, G are isomorphic. 

Indeed, for such e we can find intervals / 0 , R,..., I n such that 

^ Ii — — (®i—1) ®i+l) 

where each Jj = fl I l+ \ is nonempty . 1 Note that Jj C Jf £ C (aj,aj + i) automatically. By the 
constructibility of F, G the various inclusions of intervals induce isomorphisms 

F(/,) = F on = F(i?) = F((a j _ 1 ,a, +1 )), F(J f ) = F (J‘) = F (J?) = F ((a,,^.)), 

0(1,) = G (/') = G(/f) = G((a i _ 1 ,a i+1 )), Gf J,) = G (Jf) = G (J?) = G((o ( ,a i+1 )). 

It follows that the following maps from an e -interleaving 


F(A) 

F(i?) F(/ 2e ) 


F(Ji) 

F (Jf) 

F (J? e 


X X 

and 

X 

X 

G {Ii) 

G(i?) G(/ 2£ ) 


G (J,) ~ 

G(J?) 

= G (Jf 


1 Specifically, I; = + at) — 5, a,i + a^+i) + 5) will do, for sufficiently small S > 0. 
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induce isomorphisms F((ct;_i, a, + i)) = G((ctj_i, cq+i)) and F((aj,cti + i)) = G((a;, aj + i)). These 
isomorphisms are natural with respect to the inclusions (aj_i,aj) C (aj_i,aj + i) and (aj,a i+ i) C 
(aj_i,aj + i) because the interleaving maps are natural with respect to Jj_i C I t and Jj C I r . 

Proposition 3.12 and Corollary 3.14 imply that this collection of isomorphisms extends to an 
cosheaf isomorphism between F, G. □ 

Corollary 4.9: The Reeb distance d^ is an extended metric on isomorphism classes in Reeb. 

□ 


4.2 Smoothing functors 

We can express the notion of e-interleaving in categorical language, following [9]. Think of the 
expansion operation on intervals 


: Int —>• Int; / (->■ I £ 

as a functor (since / C J implies I £ C J £ ). Any pre-cosheaf F : Int — > Set can be ‘e-smoothed’ 
to obtain a new pre-cosheaf F12 £ : Int — > Set. Thus FO £ (J) = F(/ £ ) by definition. 

Observation 4.10: Asking for natural families of maps (</?/), {ipi ) as in equation (4.2) is pre¬ 
cisely the same as asking for natural transformations p : F G12 £ and ip : G =>- Ffi £ . 

Observation 4.11: The functor fl £ is a pointed endofunctor on Int, because the inclusions 
/ C I £ define a natural transformation lo £ : li nt 0 £ . From this we get a natural transformation 

of = Fa/ : F => FO £ 

dehned explicitly by (of); = F [I C I £ ] : F(/) —>■ F(/ £ ). 

Observation 4.12: Conditions (4.3) can be written as {ipfl £ ) ° P — and (<^0 £ ) o ip — crjf. 

The two observations combined give a more purely categorical definition of ^-interleaving. The 
restatement of (4.3) asks that the following diagrams (of natural transformations) commute: 



G 



(4.13) 


Implicitly we are using C2 £ C2 £ = 

Now we consider the smoothing operation F i —> Ffl £ in its own right. This, we claim, is a 
functor on pre-cosheaves. Indeed, any natural transformation p : F G gives rise to a natural 
transformation pfl £ : Ffl £ G0 £ dehned explicitly by {pfl £ )j = pje : F(/ £ ) —> G(J £ ). One 

verihes immediately that this procedure respects composition and identities. Thus: 

Definition 4.14 (Smoothing functor): Let S £ be the endofunctor of Pre = Set Int dehned by 
precomposition with fl £ . Thus <S £ (F) = FO £ , and <S £ [</?] = pfl £ : FhL =>■ Gfl £ for a morphism 
p : F G. 
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Observation 4.15: It follows from Observation 4.11 that each S £ is a pointed endofunctor of 
Pre, in the sense that there is a natural map dp from each pre-cosheaf F to its smoothing «S e (F). 
This is defined at each interval / to be 

F[/ C /'] : F(7) -4 F(/') = S.(F)(/). 

Succinctly, these maps comprise a natural transformation a £ : lp re =>• S e , defined of = Fca £ . 

In the remainder of this section, we study the properties of S F . The next two propositions 
support an analogy between smoothing of pre-cosheaves and smoothing in functional analysis 
(for example by convolution with a heat kernel). 

Proposition 4.16: (S £ ) is a semigroup of endofunctors (on Pre) because (f2 e ) is a semigroup 
of endofunctors (on Int): the relation Q ei+£2 = S7 £2 f2 £l implies S £1+£2 = d> ei d> £2 . P 

Proposition 4.17: Smoothing is a contraction: di(«S e (F),«S e (G)) <di(F,G). 

Proof. If (p, ib define a 5-interleaving between F, G then gjIL, AfL define an 5-interleaving between 
<S e (F) = Fn e and 5 e (G) = G Q e . D 

The next theorem indicates that we can make geometric use of S £ . 

Theorem 4.18: The functor S £ restricts to functors S £ : Csh —>■ Csh and S £ : Csh c —> Csh c . 

We split the theorem into two propositions. 

Proposition 4.19: The functor S £ carries constructive pre-cosheaves to constructible pre¬ 
cosheaves. 

Proof. Let S' be a critical set for F. We claim that S £ := (S — e) U (S + e) is a critical set for 
Ff2 e . Indeed, if / C J and J \ I does not meet S £ then J £ \ I £ does not meet S. Thus Ffl e [J C 
J] = F[/* C J £ ] is an isomorphism. And if I is contained in (—oo,min(S ,£ )) U (max(S' £ ), +oo) 
then I s is contained in (—oo, min(S')) U (max(S'), +oo) so F Q, e (I) = F(/ £ ) = 0. □ 

Proposition 4.20: The functor S £ carries cosheaves to cosheaves. 

Proof. Let F : Int — y Set be a cosheaf. We show that F1L is also a cosheaf. 

Let U be an open interval covered by open intervals (I p \ a E A). Then U £ is covered by 
(J £ | a E A). We want to show that FL2 £ (L/) is the colimit of the diagram 

]jFfi £ (J p n/ g )^]jFQ £ (J p ), (t) 

P,Q P 

knowing that F(t/ e ) is the colimit of the diagram 

II F <4 n t‘) => II F ( J h (t) 

p,g p 

The two diagrams are almost identical, except that (f) has extra terms on the left-hand side 
whenever I £ D I £ is nonempty but I p D I q is empty. We will show that these extra terms do not 
affect the colimit. 

Let Z be a set, and suppose we are given maps ( p : Ff l £ (I p ) = F(/ £ ) —>■ Z for all p, such that 

e„ o F|/« n i\ c /'] = C s o F[/' n i\ c /;] (§) 
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whenever I p D I q ^ 0. We will show that equation (§) holds in the additional cases where 
If fl If ^ 0. Then by the universal property for the colimit of (i) there will be a unique map 
C : F(£/ £ ) = F Q e (U) -)• Z such that 

C P = C° F[Ip C U £ ] = C o F n £ [i p C E7] 

and this will confirm that F £l e (U) satishes the universal property for the colimit of (f). 

To this end, let J = I p fl I £ be nonempty with I p D I q — 0, so that the two open intervals 
sandwich between them a nonempty closed interval K. Since K is compact and connected and 
contained in U, we can find a finite connected chain (I Pi ) of intervals meeting K which connects 
the two ends of K ; so I p — I po , I pi ,..., I Pn = I q with each I p . fl I Pi+1 nonempty. Now, by metric 
considerations, each thickened interval I p . contains J. Then for each i we have the following 
diagram: 


m.) 



The hve maps on the left are those assigned by F to the corresponding inclusions of intervals. 
The two triangles on the left commute since F is a functor, and the quadrilateral on the right 
commutes by (§) because I Pi fl I Pi+1 is nonempty. It follows that 

Cp, ° F [J C /'] = C K+1 O F[J C I’ t J 

so, following the chain, we deduce (§) for the pair p, q. 

In sum, we have shown that if (§) holds for all p, q with I p fl I q nonempty, then it holds for 
all p, q with I p fl I q nonempty. Thus the extra terms do not affect the colimit, and the proof is 
complete. □ 

Remark 4.21: The proof is not specific to the category Set. If F : Int —y C is a cosheaf in an 
arbitrary category C (with an initial object), then its smoothing FI2 £ is a cosheaf, by the same 
argument. 

4.3 Thickening functors 

On the geometric side there is a family of functors (%) acting in parallel to the smoothing 
functors (S e ) that act on the cosheaf side. We study these functors now. 

Definition 4.22: For £ > 0, the thickening functor % : R-Top —y R-Top is defined as follows. 

• Let (X,/) be an R-space. Then 7^(X,/) = (X e , f £ ) where X e = X x [—e, e] and f £ (x,t) = 
f(x) +t. 

• Let a : (X,/) —* (Y, g) be a morphism. Then T e [a] : (X £ ,/ £ ) —> (Y £ ,/ £ ); (x,t) (->• (a(x),t). 
It is easily conhrmed that this is a functor. 
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Observation 4.23: The thickening functor T e is a pointed endofunctor of R-Top. indeed, the 
canonical embedding of X as the zero section of X £ defines a natural transformation r £ : 1r_top =>• 
T e - Schematically we draw the picture 


and formally we define r £ = (r|) by 

Tf ■ (X, /) -t (X £ , fe)- X l-G (:X , 0). 

Naturality follows trivially from the fomula. 

The next theorem gives the precise meaning of ‘acting in parallel’: one can e-thicken before 
taking the Reeb cosheaf, or e-smooth after taking the Reeb cosheaf, and the result is the same. 

Theorem 4.24: We have CT e — S £ C. That is, the functors CT e , S £ C are naturally isomorphic. 

The main part of the proof is understanding the relationship between inverse images of / and 
f £ . Let p : X £ = X x [—£,£] —> X denote the projection onto the first factor. 

Lemma 4.25: The map p restricts to a homotopy equivalence / e ^ 1 (J) ——> / -1 (/ e ) for each 
interval I. 

Proof. Let pj denote the restriction of p to / £ _1 (J). Then pj carries / £ -1 (/) into / -1 (/ e ) because 
/(x) + f G / implies /(x) G I £ . 

To define a homotopy inverse qj : /^ 1 (/ e ) —>■ /~ 1 (/), we select a continuous function A : I s — * 
[—e,e] such that s + A(s) G / for any s G I £ . For instance, if we write / = (a, b ) then A can be 
any continuous function on (a — e, b + e) whose graph lies in the interior of the parallelogram 
shown here: 



a—s a a+e b—£ b b+£ 


There is no problem choosing such a function for any given interval I. We set qi(x) = (x, A(/(x))). 
Then qj carries /^ 1 (/ e ) into / e -1 (/), since f{x) G I £ implies f{x) + A(/(x)) G / by the condition 
on A. 

Certainly pjqj is equal to the identity on f~ l (I £ ). In the other direction, q I p I (x,t ) = 
(x,\(f(x))) and this function is homotopic to the identity on / e _1 (J) via linear interpolation in 
the f-coordinate. This works because for any fixed x the set / £ _1 (J) meets the fiber {x} x [—e, e] 
in an interval. □ 

Proof of Theorem f.2f. We will define a natural transformation p : CT e => S £ C and show that 
Pf is an isomorphism for each object / = (X, /) in IR-Top. 

Expanding the definitions, C%(f) and S £ C(f) are pre-cosheaves which evaluate on intervals 
and morphisms as follows: 

[C%(f)](I) =n 0 fp(I) [CT.(f)][I C J] = £{J lfP(I)CfP(J)) 

[■S c C(f)](I) = no /->(/*) [S«C(/)][/ C J] = lr „[/- 1 (F) C /-'(J*)] 
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We define the natural transformation pf : C%(f) =>• S £ C(f) by the formula (p/)/ = 7r o[pi\. Here 
pj is the map dehned in Lemma 4.25. Since it is a homotopy equivalence, it follows that (p/)/ is 
an isomorphism. Applying 7r 0 to the left square of the commutative diagram 


/Pm -* W) — XX [-£,£] 


PI 

PJ 



Y 



' 

-1 ( T?\ 

_^ f- 1 


_^ \ 


we confirm that pj is a natural transformation; that is, a morphism of pre-cosheaves. 

(pf)i is an isomorphism it follows that pf is an isomorphism of pre-cosheaves. 

To finish we must show that the family of pre-cosheaf isomorphisms p = (pj) 
with respect to morphisms in [R-Top. In fact, for any morphism a : (X,/) —)■ (¥, g) 
commutative diagram 

f;\i) — 9i\i) 

p { p? 

r \n ^ g~\n 

for each interval /, to which we can apply 7To to get the required naturality condition. □ 

The thickening functors % preserve constructibility. We give a simpler result first, since we 
can state its proof more briskly and it is all we need for the topological smoothing of Reeb graphs. 

Proposition 4.26: If (X,/) G Reeb then 7^(X,/) G R-Top c . 

Proof. An R-graph can be represented as a piecewise linear function / on a compact 1-dimensional 
polyhedron X. By construction, X e is a compact polyhedron and f £ is piecewise linear, so 
7^(X, /) = (X e , f £ ) is a constructible R-space. □ 

Here is the full result. 

Theorem 4.27: If (X,/) G R-Top c , then 7^(X,/) G R-Top c . 

Proof. It will be helpful to reparametrize (X e , f £ ). Consider the R-space dehned as follows: 

X £ = {(i,m)gXxR|m-£< f(x) < U + £■}, f £ (x, u ) = u. 

We claim that (X e , f £ ) is isomorphic to (X e , f E ). An inverse pair of morphisms dehned as follows: 

(Xe, fs) “>• foe, Is)', (X,t) (x,t + f(x)) 

{%, fe) (X e , / e ); (a;, u) H- (x, u - /(a;)) 

See Figure 8. For the rest of the proof, we drop the tildes and write (X e , / e ) to mean (X e , f £ ). 
Step 1: Compact locally path-connected fibers. Notice that (in the new coordinates) we have 

/rV) = r 1 ^- e,u + e\ x {«}. 

Each point of / _1 [m — e,u + e] has a neighborhood which looks like a cylinder on some E* (in 
the non-critical hbers) or a mapping cylinder to some V; (in the critical hbers). Since the Vj, E, 
are locally path-connected, the same is true for these neighborhoods. Thus each fiber is locally 
path-connected; and compact because X is compact. 


Since each 

is natural 
we have a 
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Figure 8: The thickened space X £ reparametrized as a subspace (green band) of X x IR (tall 
rectangle). The map f £ becomes the projection onto the factor IR. In this example the three 
critical values of (X, /) give rise to six critical values of (X £ , f e ). The fibers at and between the 
critical values are illustrated in Figure 9. 


Step 2: Critical set. We define S £ = {a + e,a — e | a G S}. These are precisely the values of u 
where one of the endpoints of [u — £, u + e] meets the critical set S. Away from these values, we 
find that the fibers of X £ are locally constant in topological type. Write S £ = {b 0 , bi,..., b m } in 
increasing order. 

Step 3: Critical fibers. We set VJ k = f^ 1 [b k — £, bk + e] and note that VJ k x {b k } = ffi l {bk). 

Step 4-' Non-critical fibers. We set F^ = f~ l (b) for some b G (b k ,b k +i). This fiber takes the 
following form. If S meets [b — £, b + e] in a nonempty set {dj, ctj+i,..., aj} then 

F fc = Ej_i x [b-e,ai] U / _1 [aj, afi U Ej x [a j: b + e] 

If S does not meet [b — e, b + e] then simply = E, x [b — e, b + e] for some i. 

Step 5: Cylindrical structure maps. We define : F fc x [b k ,b k+ 1 ] —>• X £ as follows. First we 
define a map a% : E k —> / _1 [u — e, u + e\ for each u G [bk, b k + 1 ]- If S meets [b — e, b + e] we use 
the following diagram: 

E k = Ej_i x \b — e, aj\ U / _1 [aj, afi U x [aj, b + e\ 

l| m I r| 

Ej_! x [u — e, a*] U f- l [a.i,aj] U E jX[aj,u + e] —> f l \ u ~ e , u + £ ] 

The map M is the identity and each map L and R is the homeomorphism defined by linearly 
stretching the second factor of the domain onto the second factor of the codomain. If S does not 
meet [b — e, b + e\ then we use the diagram 

Ffc = Ej x [b — e, b + e] ——>■ E iX[u — £,u + e] —> f~ 1 [u — £,u + £] 

where M is the homeomorphism defined by linearly translating the second factor. 

Then a k (£,u) = (ct]j(£),u) is the required cylindrical structure map. It is continuous because 
the coefficients of the stretches or translations are continuous in u. 

Note that a k is a homeomorphism when u G (b k ,b k+ 1 ). We interpret the two endpoint cases 
as attaching maps : Fand a^ fe+1 : E k —* W fc+1 . These need not be homeomorphisms. 






















C <3 > 

|oo 


Figure 9: The critical fibers VJ k (green) and non-critical fibers F fc (purple) of the thickened [R-space 
(X e , f e ) in the example of Figure 8. The fibers are homeomorphic to interlevelsets f~ 1 [a — e, a — e\. 


Step 6: Constructibility. We now have maps from the spaces W fc x {6 fc } and x [b k , b k +i\ to X £ 
which respect the attaching maps a b k k , a b k k+1 . By the quotient principle this induces a continuous 
map from the constructible [R-space built from the W&, F& and the corresponding attaching maps. 
This map is a bijection on each fiber and therefore a bijection. Since the domain is compact and 
the codomain is Hausdorff, the map is a homeomorphism. 

This completes the proof that (X £ , f £ ) is constructible. See Figure 9 for an example. □ 

Corollary 4.28: It follows from the proof that if / has critical set S then f e has critical set 
S e = (S — e) U (5 + e) = {a — e, a + e \ a e S}. 

Remark 4.29: The thickened space (X e , f £ ) can be thought of as giving a natural topology on 
the family of ‘sliding windows’ on (X, /) of width 2e. 

4.4 Topological smoothing of [R-graphs 

We are now in a position to define a semigroup (U £ ) of ‘topological smoothing’ functors on Di¬ 
graphs. We then use these functors to reinterpret the Reeb distance in a more purely geometric 
way. There are two reasonable ways to define this semigroup: 

• Use the thickening functors (77) followed by a projection onto Reeb. 

• Transfer the smoothing functors (S £ ) to Reeb using the equivalence of categories. 

We favour the first method, which gives the following explicit definition, for all e > 0: 

Definition 4.30: Define the Reeb smoothing functor U £ : Reeb —y Reeb by U £ = 1ZT £ . 

Given an R-graph / = (X, /), it follows that the fiber of U £ f over t e R can be identified with 
the set of connected components of f~ 1 [t — £,t + e\. If the vertices of the original graph occur over 
a critical set S, then the vertices of the smoothed graph occur over the set S £ = (S + e) U (S — e). 
These facts follow from Theorem 4.27 and its Corollary. 

Example 4.31: To better understand this construction, consider the examples of Figure 10. 
The initial R-graph X is given in column (a) with the function / implied by the height. In order 
to visualize the space X x [—e, e] in column (6), X is redrawn with a [—e, e] interval added at 
each point. Since these intervals are drawn vertically, the function f £ can still be visualized as 
the height function of this new space. The Reeb graph of this space is overlaid in column (c) 
and drawn by itself in column (d). 
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(1) Here the IR-graph is a line. It gets stretched by e in both directions. 

(2) The np-fork in this example gets pushed up by e. (A down-fork would get pushed down 
by e.) 

(3) This example has a loop with height (b — a) < 2e. After thickening, every levelset has only 
one connected component so the resulting Reeb graph has gotten rid of the loop entirely. 

(4) This is a more complicated mix of the ingredients above. Note that the height of the loop 
shrinks by 2e. There is interesting behavior on the right side where a down-fork interacts 
with an up-vee. 

It is easiest to access the properties of (U e ) by comparing these functors with (S £ ). 

Proposition 4.32: The functors C"U £ and S £ C" : Reeb — y Csh c are naturally isomorphic. 

Proof. We have C"1ZT £ — C'T e — S £ C" by Theorems 3.17 and 4.24. □ 

This implies, in particular, that the functors (U £ = VS £ C") suggested by the second method 
above are naturally isomorphic to the functors (W £ ). We prefer the hrst method because it is 
more geometric and because the inverse functor T> used by the second method has not been 
defined explicitly. 

Observation 4.33: The family of functors (U £ ) form a semigroup of contraction endofunctors 
(in Reeb), in the sense that: 

• W 0 - iReeb and W £1+£2 ~ U £1 U £2 for all ei,e 2 > 0; 

• d R (W £ (/),W £ (sO) < d K(f,g) for all f,ge Reeb and e > 0. 

These assertions follow immediately from the corresponding assertions (Propositions 4.16 and 4.17) 
for the family of endofunctors (<S £ ) of Csh c , thanks to Proposition 4.32. For the semigroup prop¬ 
erty, we have to replace '=’ with since that is all we can deduce from an equivalence. 

Observation 4.34: Each U £ is a pointed endofunctor of Reeb: there is a family ( £ = (Q) of 
maps 

cy.f^u £ {f) 

from each IR-graph to its e-smoothing, which constitute a natural transformation ( £ : lR ee b W £ . 
For a given graph / = (X, /), the map Q is the composite 

X -- X £ -- X £ / ~ 

where the hrst map is the inclusion of X as the zero-section of X £ and the second map is the 
Reeb quotient. That is to say, ( £ is the composite of the natural transformations r £ and pf e of 
Observations 4.23 and 2.11. Geometrically, the map sends each point of f~ 1 (t ) to the connected 
component of f~ l \t — e, t + e] that it belongs to. 

We finish this section by showing that ( £ : lR ee b =>• W £ corresponds exactly to a £ : lcsh c => <S £ . 

Proposition 4.35: For / G Reeb and e > 0, the right-hand square in the following figure 

/ C(f) C(f) 

C C[C] 

U £ (f) CK £ (f)^—S £ C(f) 
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Figure 10: Examples of smoothing. Column (a) represents the original [R-graph X with function / 
shown by vertical height. Column (6) shows the thickened space X x [—£,£] with function f £ still 
shown by the vertical height. Column (c) shows the Reeb graph of (X x [— e, e], f e ) superimposed 
on X x [—£,£]. Column (d) shows this new R-graph by itself. 
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commutes. Here C = C/ and cr = The left-hand side of the square is obtained by applying 

C to the diagram on the left, and the map at the bottom of the square is the isomorphism of 
Proposition 4.32. 

Proof. We need to verify that the square (of small functors and natural transformations) com¬ 
mutes when evaluated at an arbitrary interval /. The result of this evaluation is the left-hand 
square below. Unpacking the definitions, we find that it is the image under 7T 0 of the right-hand 
square below: 


vro f~\l) 


C[C]i 




^ f~\l) 


< 71 


VTO f~\l e ) 


f-\l) 


f7\l) 


f~\l) 


H 


R 


r\n 


The three labelled maps are inclusions. The left map L is the inclusion of / _ 1 (J) as the zero- 
section of /y 1 (/). The right map R is the inclusion of / _1 (/) as a subset of / _1 (/ e ). The 
horizontal map H is the homotopy equivalence defined in Lemma 4.25 in terms of the auxiliary 
function A. 

ft is enough to show that this right-hand square commutes up to homotopy. Indeed, this is 
the case: the maps L and HR are homotopic via a fiberwise straight-line homotopy. This may be 
understood by contemplating the figures 



which schematically represent the two homotopic maps. □ 


4.5 Interleaving of R -graphs 

We are now in a position to interpret the interleaving distance between R-graphs geometrically. 
The original version (Definitions 4.1 and 4.4) asks us to compare two cosheaves. The discussion 
in Section 4.4 allows us to do the comparison directly on the graphs themselves. Let / = (X, /) 
and g = (Y, g) be R-graphs. An e -interleaving of their Reeb cosheaves F, G is described by the 
diagrams in (4.13). We interpret this in the geometric category. 

Definition 4.36: Two R-graphs f,g are ^-interleaved if there exist maps 


a : / ~*U e g 


and /3 : g U £ f 


such that the diagrams 




and 


(4.37) 



















commute. Here a| e and /3f e are the composites 

U £ [cd ~ 

<4 : U £ f -^-*$$4 U e ^ -=-*#4<7 

UM 

PL : U e g -^-- U £ Usf -=-- U 2e f 

where the right-hand maps are the natural isomorphisms of Observation 4.33. Through the 
equivalence of categories C" and the results of the previous section, it follows that /, g are e- 
inter leaved as [R-graphs if and only if F, G are e-inter leaved as cosheaves. 

Remark 4.38: Geometrically, the map a e 2e acts as follows. Each point x E U £ f corresponds 
to a connected component of some / _1 [t — e, t + e]. The points in that component are carried 
by a to connected components of various c/ _1 [s — e, s + e] where s E [t — e, t + e]. By continuity, 
these components are all contained in a single connected component of g~ 1 [t — 2e,t + 2e\. It is 
this component that defines a 2e ( x ) E U 2e g. The map PL acts in similar fashion. 

The Reeb interleaving distance between two R-graphs is, finally, the inhmum over values of e 
for which there exists an e- interleaving between the R-graphs. 

5 Algorithms 

Our main goal in this section is to describe an algorithm for constructing the e-smoothing of 
a given Reeb graph, as well as the canonical map from the graph to its smoothing. These 
are necessary ingredients for working with the interleaving distance. This can be achieved in 
polynomial time. Calculating the interleaving distance between two Reeb graphs is not so easy: 
in general it is graph-isomorphism-hard. We can, at least, recognize interleavings in polynomial 
time. We will discuss these matters in the later subsections. 

Here is the set-up. Let / = (X,/) be an R-graph with critical set S = {a 0 , • • • , a n }. In 
computational terms, this is just a graph X with function values associated to each vertex. 
Implicitly, we assume that the function value on the edges is a strictly monotone function with 
max and min determined by the function values at the vertices. 

We wish to compute the smoothed Reeb graph U e (f). To do this, one might naively build a 
larger complex Xx [— e, e] with function f e as in column (6) of Figure 10 and then run any standard 
algorithm to compute its Reeb graph. This new complex will have one edge and two vertices 
for every original vertex, and three edges and two faces for every edge, so for a graph with m 
edges and n vertices, the new complex has 0(m + n) total simplices. Hence, the new Reeb graph 
can be computed in time 0((m + n) log(m + n)) in expectation using [27] or deterministically 
using [32], 

However, this method does not make use of the particular structure of the smoothing pro¬ 
cedure. If we do exploit that structure, we can modify the algorithm of [32] to compute U £ (f) 
while running in 0(mlog(m + n)) time. 

5.1 The smoothing algorithm 

Parsa’s algorithm for computing the Reeb graph of an arbitrary simplicial complex in [32] is a 
sweep algorithm which keeps a data structure to represent the connected components of / _1 (r) 
as the value r increases in order to determine how the connected components should be attached. 
Let / : X — » R be the original Reeb graph and f e : X x [— e, e] — >• R the thickened Reeb graph. 
Let Vi, ■ ■ ■ ,v n be the critical vertices of / sorted so that f(v i) < f{v 2 ) < ■ ■ ■ < f{v n ). To 
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simplify the explanation, we will assume general position, by which we mean f(vi) 7 ^ f(vj) and 
f[vi) ±£ 7 ^ f{vj) ±e for any i 7 ^ j. Let &i < b 2 < ■ ■ ■ b k be the sorted values {/(ty) ±e}; thus { 6 *} 
contains the critical values of f F . For the sake of notation, let [t] £ denote the interval of width 2 e 
centered at t\ that is, [ t] £ = [t — e,t + e]. Also, we will say an edge e = (v,w) with f{y) < f(w ) 
starts at v and ends at w. 

The main change here in order to instead compute the smoothed Reeb graph is the structure 
used to represent / e _ 1 (f). Essentially, using Lemma 4.25, we know that we can work equivalently 
with / j T 1 (£) or / _1 (f — e,t + e) = / _ 1 ([t] e ). Thus, rather than dealing with the larger simplicial 
complex, we determine our connected components using a graph which keeps track of edges and 
vertices within that range. We will use H for the graph data structure which represents the 
connected components of / _ 1 ([t] £ ) for the current value of t. It should be noted that there are 
three main operations required for H\ finding the particular connected component of either an 
edge or a vertex, merging two connected components, and splitting two connected components. 
This structure and methods will be throughly discussed in the next section, and the full sweep 
algorithm will be explained in the section after that. 

5.1.1 Maintenance of Level Set Representation 

As we are working combinatorially, we will consider /^ 1 ([t] £ ) as a subgraph consisting of all 
vertices with function value in [t] £ along with all edges attached to these vertices. We adopt the 
convention that we can have ‘half edges’ in this subgraph which occurs when an edge is attached 
to one vertex inside the interval and one outside. In order to minimize any confusion as we pass 
back and forth between thinking of things topologically and combinatorially, let us represent the 
topological space / _ 1 ([t] £ ) by a graph H t defined as follows: 

• every vertex in f~ 1 ([t} £ ) is represented by a vertex in Hp, 

• every edge whose interior meets /^ 1 ([t] £ ) is represented by a vertex in H t ; 

• every incidence between a vertex and an edge in / _ 1 ([t] £ ) is represented by an edge in H t . 

Notice that some edges of the original graph may only partially meet / _1 ([^] £ )- Thus, every edge 
that meets / _ 1 ([i] £ ) can have 2, 1 or 0 vertices in / _ 1 ([t] £ ). The vertices of H t that represent 
edges will therefore have degree 2 , 1 or 0 . 

Remark 5.1: The graph H t is the derived complex, or barycentric subdivision, of the part of 
the Reeb graph contained in the window / _ 1 ([t] £ ). A similar construction is implied in Figure 3. 

As in [32], we need to quickly determine the connected components of H t . This is done by 
solving the dynamic graph connectivity problem, where we store a rooted spanning forest of 
the graph H t ; this forest is the graph notated H. Since H is a subset of H t , vertices in H are 
associated to either a vertex or edge in the original graph X. 

As the value of t increases past a critical value b = f(v ) ± e, we need to be able to update 
H to continue to reflect the connected components. Algorithm 1 outlines the procedure for this, 
which will be defined as UpdateH(v, b). If b = f(v) — e, and thus b < f(v), raising the value of t 
from b — 6 to b + 5 requires adding the vertex v to H, attaching all the edges which end at v, and 
starting the edges which emanate from v. On the other hand, if b = f(v) +£, raising the value of 
t requires removing v from H, deleting the edges which end at it, and freeing the bottom of the 
edges above it. Note that vertices in H are only deleted after all attached edges are removed. 

We can now look at how to implement insert, delete, and find in this graph. In order to 
assure that we are not spending extra time adding and deleting edges, we will give each edge in 
H a weight to equal to the time that it will be deleted from H. The purpose of these weights can 
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Algorithm 1 UpdateH(u, b ) 

if b < f(v) then 

Add a vertex to H for v 
for e below v in X do 
H. insert (e, v) 
for e above v in X do 

Add a vertex to H for e 
H. insert (v, e ) 

else 

for e below v in X do 
//.delete(e, v ) 

Remove vertex e G V ( H) 
for e above v in X do 
H. delete (v, e ) 

Remove vertex v G V(H) 




Figure 11: For the portion / _1 (/) of the Reeb graph shown at left, two possible minimum 
spanning trees, H , are shown at right. In the left version, edges will only need to be removed 
from H as they are removed from / _1 (J) as I is moved up. In the far right version, edges will 
need to be added and removed as / is moved up. Thus, we prefer the example on the left to 
decrease the number of edits done to H. The edge weight in the algorithm is used to avoid the 
example on the right. 


be seen in the example of Fig. 11. A portion of the Reeb graph is at the left, and two choices for 
minimum spanning trees of the shaded region, / _1 (J), are given. Circular vertices correspond 
to vertices in / _1 (/), and square vertices correspond to edges. If the leftmost MST is chosen, 
the only edits required as the interval / moves up is to delete edges from the MST as they are 
removed from / _1 (J). If the right tree is used, the MST will require more complicated edits 
when the interval moves past the bottom vertex. Note that since every edge in H (or H t ) has 
one endpoint corresponding to a vertex and the other to an edge, we can define oj(e) to be f(y) 
for v the vertex. We will maintain that our minimum spanning tree utilizes edges with higher 
weights when possible. 

The operations needed on H (where x, aq, and aq axe all vertices in H) in order to implement 
find, insert, and delete are: 

• parent(x): return the parent of x, or null if x is the root. 

• root(x): return the root of the tree containing x. 
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• link(a;i, X 2 , w): add an edge between x\ and X 2 with weight w. 

• cut(£i,X 2 ): delete the edge between x\ and X 2 . 

• minWeight(x): return a node with minimum weight edge to its parent on the path from x 
to the root of its tree . 

• evert (x): make x the root of its tree. 

The dynamic graph connectivity problem is well studied with methods that include RC-trees [1] , 
link-cut trees [38,39], top trees [3,42], and sparsification [22], These can be implemented in order 
to perform the above operations in 0(log n) worst case, amortized, or expected time where n is 
the number of vertices in the graph; thus we leave the extended discussion of these specifics to 
the references. 

Given the above methods, we can implement the three major operations as follows: 

• find(x): 

return root (a;) 

• insert (e): 

w = oj(e) for edge e = x±x 2 
if root(xi) =root(a: 2 ) then 
evert (xi) 

x',w' =minWeight (x 2 ) 
if w' < w then 

cut (x', parent (x')) 
link(xi, X 2 , w) 

else 

link(a:i, x 2 , w) 

• delete(e): 

if e E E(H) then 
cut(e) 

Notice that since the first set of operations can be implemented in O(logn) time, find, insert, 
and delete can be as well. 

5.1.2 Full Algorithm 


Algorithm 2 Sweep algorithm 
Set H to be an empty graph 
for i — 1, • • • , k where bi = f(vj ) ± e do 
L c = Lower Comps (vj, bi) 

UpdateH(u, b) 

U c = UpperCornpsfTj, b,) 

UpdateReebGraph(L c , U c ) 


The pseudocode for the sweep algorithm is given in Alg. 2. Here, we work our way up 
the potential critical values bi, keeping track of the change in H t and using this to build the 
smoothed graph g : Y —» IR. For any noncritical t, the components of H t are associated to an 
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edge in Y whose lower vertex v and upper vertex w satisfy g(v) < t < g{w). We will keep track 
of this association by pointing the representative of a component in H to the lower vertex of its 
associated edge in Y. 

At the beginning of a step, H represents the connected components for H t with t = — S 

for a sufficiently small 5. First, we fold the components in H that could be impacted by the 
addition or deletion of Vj and associated edges using the LowerComps (vj , b {) subprocess, Alg. 3. 
These are stored in L c . Note that these are exactly the edges which end at u, where v is the 
vertex added at function value 6* in the smoothed Reeb graph. 


Algorithm 3 LowerComps (u, b) 


if b < f(v ) then 

for edges e ending at v do 
c = if.fold (e) 
if c is not marked then 
L c .add(c) 

Mark c as listed 


else 


L c = if.foid(u) 


Then the H graph is updated so that it now represents the connected components for t = bi+5 
using the UpdateH(u,&) subprocess, Alg. 1 discussed in the previous section. The components in 
the new H are determined using UpperComps(uj, 5;) (which is symmetric to LowerComps^, bi) 
and therefore not repeated here) and stored as U c . These components are the edges which start 
at v. 

Finally, to update Y, we use UpdateReebGraph (Alg. 4) to add a new vertex v to the graph. 
An edge is added for each lower component in L c which starts at the associated start vertex and 
ends at v. Then each component in U c is assigned v as the start vertex. 


Algorithm 4 UpdateReebGraph(L c , U c ) 

if #|L C | = #|C/ c | = 1 then 
return 
else 

Create a new node v in Reeb graph 

Assign that node to all c EU C 

Add an edge between v and v c for all c e L c 


5.2 Analysis of the smoothing algorithm 

We show that the overall running time is linear in the total number of simplices of the origi¬ 
nal Reeb graph. In particular, for a graph with n vertices and m edges, the running time is 
0(mlog(n + m)). First, note that the number of vertices in H is at most n + m, thus the time 
for any find, insert, or delete is 0(log(?n + n)). Every edge of the original graph is used once 
for LowerComps and once for UpperComps. Since each of these utilizes one fold operation, the 
total time spent in them is 0{m\og{n + m)). Likewise, edges are added to H twice for each 
original edge, and these edges are each deleted. Thus the total time spent in UpdateH is also 
0(m log(n + m)). Thus, the smoothing algorithm runs in 0(m log(n + m)) time. 
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5.3 Morphisms between graphs with different critical sets 

In recording morphisms, it is helpful to notice that we can get away with storing less information 
than is indicated by Proposition 2.8. In particular we can to write down maps between Reeb 
graphs without needing to subdivide the graphs to have a common set of critical values. 

Consider a map y : (X,/) —> (Y, g) given by data ipf ■ Lf —> V Y and ipf : Ef —)■ Ef as 
specihed in in Proposition 2.8. Say (Y, g) is a subdivision of (¥', g). Thus |Y| = |Y'| but Y' may 
have some vertices with up and down degree both 1 which are not in Y. We abuse notation by 
calling both maps g. Since we assume that g is monotone on any edge, if there is a vertex in X 
which maps to v G V Y under y, the same information is stored if we say that it maps to the edge 
e G Ef . Thus, for the non-subclivided version of </?, we instead have a map tpf : V) x —> V Y Li Ej . 
In order to get rid of confusion with indices, we will denote this map 

tpV :V*^V Y UE Y ' 

while remembering that a vertex v G X will map to either a vertex w G V Y ' with f(v) = g(w) or 
to an edge e = ( u,w ) G E Y ' with g(u ) < f(v ) < g(w). 

Likewise, the image of an edge e = ( u , v) G Ef in Y' is a monotone path which begins at 
ip v (u) and ends at (p v {y). In terms of the combinatorial structure, this path is a sequence of 
edges ei, e 2 , ■ ■ ■ , e*, such that the top vertex of e* is the bottom vertex of ej. 

5.4 The canonical map from an R-graph to its smoothing 

Let (X,/) be an [R-graph and let (Y, g) = (X £ , f £ ) be its smoothing. In order to recognize 
interleavings we need access to the canonical map ( = Q : (X,/) —> (Y ,g). We show how to 
obtain this in terms of the alternate description of maps described in the preceding subsection. 

We do this during the course of constructing (Y, g). During the sweep, in addition to stopping 
at critical values {6*} of the function g we will also stop at the critical values {a*} of /. This way, 
we can determine £( v ) by determining the component of E[ containing v, and finding the edge 
or vertex of Y to which it is associated. (It is usually an edge, unless by chance we have some 
Icq — cij | = e.) We retain this information for each vertex v G R(X). For each edge e G -E(X), 
we must record the corresponding path C( e )- This we do by maintaining a list with the edge. 
At each update of //, we find the representative of the component of E[ containing e G V(H). 
Thus, when the algorithm ends, the map ( has been completely determined. 

5.5 Complexity of the Reeb interleaving distance 

We have not given a method for calculating the Reeb graph interleaving distance, because it is 
not easy. Let us first take the inhmum out of consideration, by selecting e > 0. This leads to 
the question in the next result. 

Proposition 5.2: “Can (X,/) and (Y, g) be £-interleaved?” is in NP. 

Proof. We need to identify possible certificates and a polynomial-time verification that those 
certificates guarantee an e-inter leaving. We begin by constructing U e f , W 2e /, U £ g, U 2e g and the 
morphisms and (ff. A candidate certificate is a pair of maps «:/—>• U £ g and f3 : g —> U £ f. 
We calculate a\ £ and /3f e and the composites a\ £ o /3 and /3f e o a and return the answer “yes” if 
the equations 

<4 ° P = Cf and /4 o a = (f 

both hold. By (4.37), the correct answer is “yes” if and only if there exists such a pair a,/3 
satisfying these tests. All the constructions and checks can be done in polynomial time, so the 
problem is in NP. □ 
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This theorem leaves open-ended the difficulties of hireling such an interleaving. When e = 0 
we have the following obstruction. 

Proposition 5.3: “Is there a O-interleaving between (X,/) and (Y, < 7 )?” is graph-isomorphism 
hard. 

We thank Tamal Dey and Jeff Erickson for the following argument. 

Proof. We find a reduction from the graph isomorphism problem to the O-interleaving problem. 
Let G and H be two finite graphs for which we wish to test isomorphism. Let us vertices in V(G) 
and V ( H) as basepoints. Is there an isomorphism G = H which sends basepoint to basepoint? 

To answer this question, convert each graph to an IR-graph by using the distance from base- 
point as a function. Some edges may need to be split in two (if their two vertices are equidistant 
from the basepoint, meaning that the distance function will increase towards the middle of 
the edge); otherwise the vertices and edges remain the same. Then these two IR-graphs are 0- 
interleaved if and only if they are isomorphic as IR-graphs (by Proposition 4.7), if and only if 
G = H preserving basepoints. 

To test whether G, H are isomorphic without reference to basepoints, it suffices to repeat this 
test keeping the basepoint of G fixed and varying the basepoint of H over all possible vertices. 
If the test fails every time, then G , H are not isomorphic to each other. If it succeeds even once, 
then they are isomorphic. I 11 this way, a solution to the O-interleaving problem for IR-graphs gives 
a solution to the graph-isomorphism problem. □ 

By Proposition 4.8, it follows equivalently that “Is the interleaving distance between (X, /) 
and (Y, g) equal to zero?” is graph-isomorphism hard. 

6 Discussion 

There is a sense in which the development in this paper is self-annihilating: all of our cosheaf 
constructions are in the end realized geometrically. At least, this is true in the constructible case. 
With that in mind, there are perhaps two main motivations for our work: 

• The correspondence between constructible cosheaves and IR-graphs allows us to transfer 
ideas from one realm to the other. For instance, the Reeb cosheaf distance occurs very 
naturally in the context of persistence but its geometric equivalent is not an obvious con¬ 
struction. 

• The realm of cosheaves is broader that the realm of IR-graphs, since we can work quite 
easily with non-constructible cosheaves whereas IR-graphs are necessarily constructible. We 
don’t claim any immediate applications for this greater generality, but it is good to know 
that it is available. 

A pleasant consequence of this thinking is that the Reeb graph emerges as yet another instance 
of topological persistence, standing alongside the persistence diagram and the dendrogram (or 
join-tree) as a persistent invariant. It shares with those invariants an ‘interleaving’ strategy for 
defining a distance, and an easily accessed stability theorem. 

Here are a few closing remarks. 
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Persistence. Our approach was inspired by the paper of Morozov, Beketayev and Weber [30] 
who defined interleavings of join-trees (called ‘merge trees’ in that paper) in a geometric way. 
This ties in nicely with Bubenik and Scott’s approach to persistence [9]: their ‘generalized 
persistence modules’ are functors from the real line to an arbitrary category, and interleavings 
can be defined in terms of translations of the real line. By regarding a join-tree as a Set-valued 
persistence module, the two approaches lead to the same interleaving distance and the same 
stability theorem. 

More generally, one may define interleaving distances for persistence modules over an arbitrary 
poset [ 8 ] . In our present work, we regard Reeb graphs as Set-valued functors on the poset Int of 
open intervals in the real line. Because these functors are cosheaves, we get a tight relationship 
between the geometric and the category theoretic points of view [24,43,48]. That said, the 
cosheaf condition is not at all needed to define the interleaving distance or to obtain the stability 
theorem. Functors on Int analogous to the Reeb cosheaf may be obtained by replacing 7 r 0 with, 
for example, 717 or hR. Since the result is usually not a cosheaf, the question remains how to 
manage these objects. 

Higher dimensions. Reeb graphs easily generalize to Reeb spaces. Given a continuous map 
/ : X — > M, we say two points i,n G X are equivalent if f(x) = f(y ) = p and if x and y lie on 
the same path-component of the fiber / _ 1 (p). The Reeb space is the resulting quotient space 
together with its induced map to IM. 

The structure of constructible Reeb spaces is rich and much is known. For example, suppose IM 
is a manifold and suppose the map / is stratihable with respect to a stratification S of IM. Then its 
Reeb space is also stratihable with respect to S. Using an argument similar to the one presented 
in this paper, one can show that an iS-constructiblc Reeb space is equivalent to a hnite-set-valued 
5-constructible cosheaf over IM. 

Going further, Robert MacPherson observed (in unpublished work) that 5-constructible 
cosheaves can be completely classified in the following way. He constructed a category Ent(IM, S ), 
called the entrance path category , whose objects are the points of IM and whose morphisms p —> q 
are homotopy classes of entrance paths from p to q. These are paths 7 : [0,1] —> IM from p 
to q which are allowed to move only from higher- to lower-dimensional strata. More precisely, 
whenever s < t, the stratum containing 7 (t) is contained in in the closure of the stratum con¬ 
taining 7(s). Homotopies are required to remain within this class of paths, with fixed endpoints. 
MacPherson’s result is that 5-constructible cosheaves are essentially the same thing as func¬ 
tors on Ent(IM,5). A special case is seen in Propositions 3.18 and 3.12, where constructible 
cosheaves with critical set S and their morphisms are shown to be equivalent to functors on the 
zigzag diagram (3.10) and natural transformations between them. Here, diagram (3.10) is inter¬ 
preted as a category, and this category is equivalent to Ent(IR, S ). We recommend the papers of 
Treumann [43], Woolf [48] and Curry [17] as further reading on these topics. 

Edelsbrnnner, Harer and Patel [21] gave the first algorithm to compute the Reeb space of 
a piecewise linear map from a simplicial complex to the plane. Their algorithm computes the 
coarsest stratification of the Reeb space. There is much room for improvement in its running 
time. 

Non-Constructibility. In onr treatment, we assign special importance to the constructible 
objects in the categories of R-spaces and cosheaves. We do this partly for the sake of the 
equivalence theorem: we do not know how to interpret non-constructible cosheaves in a cleanly 
geometric way. And the loss of generality is not too great, because many regular situations (Morse 
functions on a compact manifold, definable functions on a compact semialgebraic set, etc) supply 
us with constructible R-spaces. Our definition is perhaps too broad: our constructible R-spaces 
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are built out of spaces which are locally path-connected, which is enough for our purposes. On 
the other hand, if we want to use higher-dimensional functors n *. or we ought to request local 
contractibility. The usual well-behaved examples have this stronger property. 

On the other hand, much of the theory works just as well for non-constructible cosheaves; so 
it would be a pity to rule them out altogether. Is there a corresponding equivalence theorem? 
One promising possibility is to define Stone-valued rather than Set-valued Reeb cosheaves. The 
connected components of a topological space can be regarded not merely as a set, but as a 
quotient space. If the original space is compact, then the quotient is both compact and totally 
disconnected, and is called a Stone space. Onr suggestion is to consider R-spaces (X, /) where X 
is locally compact and Hausdorff, and / is proper. Then is a Stone space for every 

compact interval [a,b\. In this setting, Example 3.20 is no longer a point of failure of the cosheaf 
condition for ttq: the colimit described there becomes a non-Hausdorff two-point space when 
evaluated in the topological category, and therefore becomes the desired one-point space when 
evaluated in a category of Hausdorff spaces such as Stone. All of this can be done in the category 
of (not necessarily compact) totally disconnected spaces, but we suspect that Stone spaces may 
be more fruitful; in part because their category is known to be dual to the category of Boolean 
algebras. A necessary requirement is to work with cosheaves over the ‘site’ of compact intervals 
rather than the more usual open intervals. 

Simplification versus smoothing. We have deliberately chosen not to describe our topolog¬ 
ical smoothing algorithm for Reeb graphs as topological ‘simplification’, a term that is widely 
used elsewhere [34]. The distinction between simplification and smoothing is as follows. Simpli¬ 
fication seeks to remove unnecessary loops and edges; and the methods that exist carry this out 
in an ad hoc local fashion by, for instance, collapsing small loops. Larger features are expected to 
retain their metric properties; this makes simplification a difficult problem. Smoothing, on the 
other hand, modifies the graph and its function on a global scale. There are simplifying effects, 
such as the disappearance of small ‘noisy’ loops as the smoothing parameter is increased; but 
there are also global side-effects such as the steady divergence of the maximum and minimum 
values of the function. We haven’t in this work attempted to identify what specific geometric 
simplifications take place as a result of smoothing. Some such questions are addressed in [5]. 

Computational complexity. While the computation of the smoothed Reeb graph is algorith¬ 
mically reasonable, we are still working to find sensible strategies for computing or estimating the 
interleaving distance. We can narrow our concerns to specific values of e, using binary splitting 
to converge towards the true value. For a fixed e, the existence of an e -interleaving we have seen 
to be difficult. The obvious naive algorithm would follow Proposition 5.2, testing all possible 
morphism pairs a : / —> U e g and f3 : g —> U £ f for their validity as certificates. This will certainly 
have worst-case exponential running time. 

That said, we think there are reasons to be hopeful. We know that e = 0 is difficult, but on 
the other hand when e = large all the internal structure disappears and the question becomes 
trivial. What happens along the way? There may be classes of R-graphs which can be compared 
efficiently; and there may be worst-case bad algorithms which behave well in practice most of 
the time. We expect that there exist easily-computed invariants that obstruct the existence of 
e- interleavings; a trivial example is given in Proposition 4.7, and there should be many more. 
There may be ways of selecting certificates a, /3 intelligently, taking the required relations into 
account. It seems to us that there are plenty of open questions here. 

It has also been shown [5] that the interleaving distance is tightly related to another metric 
on Reeb graphs, the functional distortion distance [4]; in particular, the two metrics are strongly 
equivalent. It is possible that this relationship will allow for an interplay of methods for com- 
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putation between the two methodologies; however, it is more likely that this relationship will be 
most directly useful for passing any NP-hardness results between the two metrics. 

We finish by acknowledging our other goal in writing this paper: to support the introduction 
of category theory and sheaves to the computational geometry community. There have been 
many other such efforts which we are proud to stand alongside. To name a few: Robert Ghrist, 
Justin Curry, Sanjeevi Krishnan, Michael Robinson and Vidit Nanda have developed many new 
applications of sheaves, cosheaves and their generalizations; and Peter Bubenik and Jonathan 
Scott have introduced category-theoretic thinking to the study of persistence. We hope that our 
direct presentation of cosheaf methods will help promote their work. Sometimes a Reeb graph 
is just a Reeb graph; but very often, secretly, there is a cosheaf waiting to get out. 
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